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ABSTRACT 
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INTRODUCTION 


I.  SOLID  STATE  DEVICE  RESEARCH 

Multi-energy  Be  implants  in  n-type  InP  and  In.  Ga  As  P,  have  yielded  layers 

1 "x  x y 1 ~y  ao  _ -2 

with  uniform  as-implanted  Be  concentrations  of  approximately  3 x 10  cm  . 

14  -2 

Sheet  carrier  concentrations,  which  were  in  the  (1-2)  x 10  cm  range,  were 
obtained  on  samples  implanted  at  room  temperature  and  annealed  at  750°  and 
700 °C  for  InP  and  InGaAsP,  respectively. 

P-n  junction  diodes  have  been  formed  by  Be- implantation  in  n-type  InP  and 

In,  Ga  As  P,  . Low  leakage  currents  and  abrupt  voltage  breakdowns  were 
i—x  x y i*y 

observed  for  both  mesa  and  planar  InP  diodes.  Scanning  photoresponse  measure- 
ments of  the  In,  Ga  As  P,  (E  * 1.0  eV)  mesa  diodes  showed  uniform  ava- 
i-x  x y l-y  g 

lanche  gains  of  2 to  3 times. 

Liquid-phase  epitaxy  has  been  used  to  reproducibly  grow  InP  and  InGaAsP  alloys 

15-3 

with  Np  - NA  at  the  low  10  -cm  level.  The  77  K electron  mobilities  for  the 
InP  are  in  the  40,000-  to  60,000 -cm2/v-sec  range,  with  Njj/N^ between  2.5  and  6. 
For  tag  gGag  2^so  5P0  5'  77  ^ mobilities  are  12,000  to  14,000  cm2/V-sec 

with  ND/NA^  2. 

Analysis  of  measurements  of  the  high-magnetic-field  Hall  coefficient  vs  tempera- 
ture yields  values  for  the  concentrations  of  deep  donors  in  InP  and  InGaAsP 
specimens  of  5 x 10*^  cm  ^ and  3 x 10*^  cm  respectively.  The  transport  data 
were  fit  using  deep  level  donor  activation  energies  of  0.29  and  0.12  eV  for  InP  and 
InGaAsP,  respectively.  These  energy  values  were  inferred  from  the  photolumi- 
nescence spectra. 

In  a study  of  photoluminescence  of  LPE-grown  InP  and  InGaAsP,  several  spectral 
peaks  were  observed  in  a range  0.3  to  1.41  eV.  An  investigation  is  presently 
under  way  to  test  the  hypothesis  that  one  of  the  peaks  in  the  photoluminescence 
spectra,  with  peak  energies  of  1.12  eV  for  InP  and  0.9  eV  for  InGaAsP,  is  due  to 
oxygen. 

The  oscillatory  interband  magneto-transmission  has  been  measured  on  samples 

of  LPE-grown  InQ  77GaQ  23^-so  52P0  48  layers-  The  band  parameters  determined 

from  the  analysis  of  the  data  are  E = 1.07  eV,  m * = 0.061  m , and  E = 17.6  eV, 
J g c o’  p ' 

where  these  parameters  are  the  direct  energy  gap  at  T * 20  K,  the  conduction- 
band  effective  mass,  and  the  k • p interaction  energy,  respectively. 


II.  QUANTUM  ELECTRONICS 

4 4 

The  F3/2  — Ifj/2  emission  cross  section  of  NdP,.014  at  1.32  pm  has  been  de- 
termined from  both  spontaneous  and  stimulated  emission  measurements.  The 
relative  performance  of  a flashlamp-excited  NdP5014  laser  at  1.32  and  1.05  pm 
has  also  been  evaluated.  At  1.32  pm,  a multimode  output  energy  of  2.7  mJ  has 
been  achieved  with  1.6  J input. 
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A CW-pumped  NLMgF-,  laser  operating  at  1.62  pm  has  been  repetitively 
Q-s  witched  at  100  Hz,  yielding  a peak  power  output  of  20  W with  a power  enhance- 

-a  .3  -4 

ment  of  >10  . A distributed  loss  coefficient  $ 10  cm  has  been  measured  in 
the  NiiMgF^  laser  crystal. 

Phasematched  fourth-harmonic  generation  in  CdGeAs,  has  been  carried  out.  The 
measured  effective  fourth-order  nonlinear  coefficient  is  X,  Jff/c  = 8 x 

"27  3 ^ i#cii  o 

10  m / v . The  contribution  of  cascade  processes,  connected  with  lower- 
order  nonlinearities,  to  the  effective  third-  and  fourth-order  nonlinear  coeffi- 
cients has  been  considered. 

A tunable  sideband  submillimeter  spectrometer  has  been  constructed  as  a first 
step  in  developing  submillimeter  frequency  standards.  Several  rotational  transi- 
tions in  DzO  and  CO  have  been  measured  to  high  accuracy, 

III.  MATERIALS  RESEARCH 

In  order  to  investigate  the  usefulness  of  laser  heating  for  removing  ion- 
implantation  damage  in  GaAs,  a study  has  been  made  of  the  effects  of  annealing 
with  a CW  Nd:YAG  laser  on  the  electrical  properties  of  Se-implanted  samples. 
Good  electrical  activation  of  the  Se  donors  has  been  achieved,  but  thermal  stress 
due  to  laser  heating  can  result  in  the  formation  of  (111)  slip  planes  that  degrade 
the  electrical  characteristics. 

To  determine  why  Cr  doping  of  InP  has  not  yielded  resistivities  as  high  as  those 
obtained  for  Fe-doped  semi-insulating  crystals  used  to  provide  substrates  for 
epitaxial  growth,  the  activation  energy  for  thermally  exciting  electrons  from  the 
Cr  acceptor  level  to  the  conduction  band  of  InP  has  been  determined  from  mea- 
surements of  the  Hall  coefficient  as  a function  of  temperature.  The  activation 
energy  was  found  to  be  0.39  eV,  indicating  that  Cr-doped  samples  have  lower  re- 
sistivities because  the  Cr  level  is  considerably  closer  to  the  conduction  band  than 
is  the  Fe  level  (which  has  an  activation  energy  of  0.65  eV),  not  because  Cr  is  less 
soluble  in  InP  than  Fe. 

IV.  MICROELECTRONICS 

A GaAs  monolithic  integrated  circuit  capable  of  receiving  submillimeter- 
wavelength  radiation  coupled  through  a high- resistivity  GaAs  substrate  has  been 
fabricated.  The  circuit  consists  of  a slot  antenna  coupled  to  a surface-oriented 
mixer  diode  by  a section  of  coplanar  transmission  line  and  an  integrated  bypass 
capacitor.  The  conversion  loss  measured  with  this  mixer  module  at  350  GHz 
was  approximately  20  dB. 

A system  has  been  designed  to  make  noise  measurements  related  to  charge  trans- 
fer and  charge  detection  on  the  100-  x 400-element  imager  for  the  GEODSS  pro- 
gram. Using  double-correlated  sampling,  noise- equivalent  signals  of  10  electrons 
at  a 400-kHz  data  rate  have  been  achieved  from  the  floating  diffusion  output  cir- 
cuit of  the  imager. 
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b: 


A theoretical  analysis  has  been  made  of  the  effect  of  sampling-finger  width  on  the 
frequency  response  of  the  SAW/CCD  buffer  memory,  and  this  analysis  has  been 
compared  with  the  measured  response  of  a prototype  device.  The  measured  re- 
sponse of  the  device  falls  off  as  predicted  by  theory. 

Improved  process  control  for  parallel-plate  plasma  etching  of  low-pressure 
chemical  vapor  deposited  films  of  silicon  nitride  and  polycrystalline  silicon  has 
been  achieved  by  constraining  wafer  temperatures  to  near  room  temperature. 
The  lower  etching  rates  obtained  with  this  method  are  not  significant  when  com- 
pared with  the  improved  uniformity  and  reproducibility  over  methods  that  allow 
the  wafer  to  be  heated  by  the  plasma. 

V.  SURFACE-WAVE  TECHNOLOGY 

A simplified  model  for  the  charge-storage  process  in  acoustoelectric  surface- 
acoustic-wave  devices  has  been  adopted;  the  many  types  of  parametric  interac- 
tions which  occur  in  these  devices  can  then  be  summarized  within  a unifying  con- 
cept. Further,  this  conceptual  view  also  identifies  the  kinds  of  signal- processing 
functions  which  are  attainable.  Devices  which  employ  two- signal  parametric  in- 
teractions are  straightforwardly  interrelated,  and  the  concept  also  extends  to 
higher-order  interactions. 

A series  of  very  sharp  stop  bands  has  been  observed  in  the  transmission  response 
of  a 200-groove  normal- incidence  grating  when  the  wavelength  of  an  incident  sur- 
face acoustic  wave  is  within  ±20  percent  of  the  grating  period.  These  stop  bands 
are  produced  by  resonant  conversion  in  the  grating  of  surface  waves  into  bulk 
waves  which  reflect  back  and  forth  from  bottom  to  top  of  the  crystal. 
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I.  SOLID  STATE  DEVICE  RESEARCH 


A.  BERYLLIUM-ION  IMPLANTATION  IN  InP  AND  In,  Ga  As  P, 

1 —x  x y l — y 

Ion  implantation  will  almost  certainly  play  a significant  role  in  the  fabrication  of 

optoelectronic  devices  in  InP  and  In,  _ Ga  As  P,  . Some  preliminary  results  on  both  n-  and 

x y i -y  1-4 

p-type  ion-implanted  layers  in  InP  have  already  been  reported.  In  this  section  we  will  pre- 
sent more-detailed  information  on  Be-ion-implanted  InP  and  report  for  the  first  time  ion- 
implantation  results  in  In^_xGaxASyPj_y.  The  results  in  InP  include  the  effects  of  implant 
and  anneal  temperature  on  the  electrical  characteristics  of  the  implanted  layer  and  the  effects 
of  implanted  Be -concentration  on  junction  depth  and,  therefore,  presumably  indiffusion  of  the 
implanted  Be.  We  will  also  report  results  on  the  use  of  a multi-energy  Be -implant  schedule 
designed  to  create  a uniform -concentration  p-type  layer  in  both  InP  and  Int  ^Ga^s^Pj  . 

The  latter  results  are  applicable  to  the  fabrication  of  high-quality  p-n  junctions,  which  are  re- 
ported in  Sec.  B below. 

The  InP  samples  used  in  these  experiments  were  cut  from  pulled  single  crystals.  High- 

resistivity  (111)  Fe-doped  samples  were  used  for  initial  experiments,  and  (100)  n-type  samples 
16  —3 

(n  « 2 x 10  cm  ) were  used  for  p-n  junction  measurements.  After  cutting  and  polishing,  the 

S 

InP  samples  were  etched  in  1:1;5;1  mixture  of  HAc:HC10, :HNCK:HC1.  The  In,  Ga  As  P, 

^ J 1 -x  x y i -y 

samples  used  were  n-type  [n  « (0.5  to  1)  x 10  cm  ] epitaxial  layers  grown  by  liquid-phase 
epitaxy  on  (100)  n+-InP  substrates.^ 

Beryllium  was  implanted  into  the  etched  surfaces  of  the  InP  samples  and  into  the  as -grown 

surfaces  on  the  In,  Ga  As  P,  samples.  Following  implantation,  two  slightly  different 
i -X  x y i -y 

methods  of  encapsulation  and  anneal,  which  gave  identical  results  on  both  implanted  and  control 
samples,  were  employed.  Both  methods  rely  on  a pyrolytic  phosphosilicate  glass  (PSG)  encap- 
sulation, nominally  7.84  wt%  P,  deposited  at  2 50°C.  Details  of  these  methods  and  first-anneal 
sample  processing  can  be  found  in  Ref.  7. 

For  the  first  set  of  experiments,  400-keV  Be+  ions  were  implanted  into  Fe-doped  sub- 
strates at  temperatures  ranging  from  room  temperature  to  200°C.  Following  implantation, 
the  anneal  temperatures  used  ranged  from  500°  to  750°C.  Figure  1-1  shows  the  sheet  hole 

concentration  and  mobility  vs  anneal  temperature  for  samples  implanted  at  150°C  with 
14-2  + 

1X10  cm  400-keV  Be  ions.  Data  for  several  samples  implanted  at  room  temperature 
are  also  shown  and  are  in  substantial  agreement  with  150°C  data.  As  shown,  there  is  a sig- 
nificant increase  in  the  Be+  activation  between  500°  and  600°C.  Between  600°  and  700°C,  the 

sheet  hole  concentration  increases  from  2.6  x 1013  to  5.3  x 1013  cm-2.  Following  a 750°C 

13  -2 

anneal,  the  sheet  concentration  in  samples  implanted  at  150°C  was  5.6  x 10  cm  , while  in 

13  -2 

samples  implanted  at  room  temperature  it  was  slightly  higher  at  5.9  x 10  cm  . 

Figure  1-2  shows  the  sheet  hole  concentration  and  mobility  vs  implant  temperature  for 
samples  implanted  with  several  different  doses  of  400-keV  Be+  ions  and  annealed  at  750°C. 

The  samples  implanted  with  3 x 1013  and  1 x 1014  cm"2  showed  only  a slight  reduction  in  con- 

14  -2 

centration  with  implant  temperature.  With  a 3 x 10  cm  implant,  however,  the  sheet  con- 
centration in  samples  implanted  at  room  temperature  was  about  a factor-of-2  lower  than  in 
samples  implanted  at  >100°C. 
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Fig.  1-1.  Sheet  hole  concentration  and  mobility  Fig.  1-2.  Sheet  hole  concentration  and  mobility 
vs  anneal  temperature  for  InP  implanted  with  vs  implant  temperature  for  InP  implanted  with 


i x lO1^  cm  2 400-keV  Be+  ions 


400-keV  Be+  ions. 


LSS 

JUNCTION 
DEPTH  \ 


ROOM 

TEMPERATURE 

9 


* / 

/ / 
o 


' 

r /a 

i 


I / Be  IN  In  P 

I / IMPLANT  SINGLE-ENERGY  MULTI-ENERGY 

o a TEMPERATURE  CC)  400  keV  FLAT  PROFILE 


2.0  SO 

JUNCTION  DEPTH  (pm) 


Fig.  1-3.  Junction  depth  vs  maximum  implanted  Be  concentration  for  400-keV 
single-  and  multi-energy  Be  implants  in  InP  and  In,  Ga  As  P,  . 


To  determine  if  these  differences  could  be  due  to  a Be-concentration  dependent  diffusion 

effect  as  has  been  observed  in  GaAs  (see  Refs.  8 through  11),  the  junction  depths  in  a series  of 

n-type  InP  |n*2x  lO1^  cm"J)  samples  implanted  with  both  a single-energy  (400-keV)  and  a 

multi-energy  Be  implant  schedule  were  determined.  The  multi-energy  schedule  was  used 

since,  in  the  actual  fabrication  of  p-n  junction  diodes,  some  type  of  multiple  implant  would 

most  likely  be  used  to  avoid  the  possibility  of  a buried  p-type  layer.  The  multi-energy  sched- 
12  13 

ule  ’ used  was  designed  to  give  a fairly  uniform  as-implanted  Be  concentration  and  was 

based  on  N at  400  keV,  12/18  N at  200  keV,  7/18  N at  100  keV,  and  4/18  N at  50  keV,  where 

14  -2 

N is  the  dose  at  400  keV.  For  N = 1.8  X 10  cm  , the  uniform  as-implanted  Be  concentration 
18  -3 

is  about  3 x 10  cm  . Following  anneal  and  PSG  removal,  the  samples  were  coated  with 

O 

wl000  A of  SiC>2  and  then  cleaved  and  preferentially  etched  in  a KOH,  KjFefCN)^  (potassium 
ferricyanide)  solution  to  reveal  the  position  of  the  p-n  junction.  The  SiC^-InP  interface  was 
used  as  the  surface  reference.  The  error  in  junction  determination  is  estimated  to  be  about 
0.2  |un. 

Figure  1-3  shows  the  junction  depth  vs  maximum  implanted  Be  concentration  (Np)  in  InP 

samples  implanted  with  either  a single-  or  multi-energy  implant  schedule  and  annealed  at  750°C. 

The  maximum  implanted  Be  concentration  was  calculated  from  the  implanted  dose  using  LSS 
12  13 

range  theory.  ’ Also  shown  is  the  junction  depth  in  several  In.  Ga  As  P,  samples  im- 

i “X  x y l — y 

planted  with  a multi-energy  implant  schedule  and  annealed  at  700°C.  Note  that  one  of  the  InP 

data  points  falls  right  on  top  of  an  In.  Ga  As  P,  data  point.  The  heavy  dashed  line 

i-x  x y l-y  ^2  33  36  .3 

represents  the  junction  depth  predicted  from  LSS  theory  ’ in  n *2  x 10  cm  InP.  The 
predicted  junction  depths  for  the  quaternary  samples  are  very  close  to  this  line  (within  0.05  pm). 
The  light  dashed  lines  are  drawn  through  the  single-energy  room -temperature  and  150°C  im- 
plant results.  Except  for  possibly  the  lowest  dose,  3 x 10**  cm”2  (Np  = 4.7  x 10*^  cm"*),  the 
junctions  depths  in  samples  implanted  at  100°  and  150°C  are  significantly  deeper  than  in  similar 

samples  implanted  at  room  temperature.  For  the  samples  implanted  at  room  temperature  with 

14  -2  18  -3 

a single-energy  implant  dose  41  x 10  cm  (Np  » 1.6  x 10  cm  ),  the  measured  junction 
depths  were  identical  to  those  predicted  by  LSS  range  theory,  indicating  insignificant  indiffu- 
sion of  the  implanted  Be.  Samples  implanted  at  room  temperature  with  a multi-energy  sched- 

18  —3 

ule,  which  resulted  in  an  as-implanted  concentration  of  3 x 10  cm  , exhibited  junction  depths 

which  were  within  experimental  error  of  that  predicted  by  LSS  theory.  In,  Ga  As  P,  sam- 

1 “X  x y » -y 

pies  implanted  to  this  level  also  showed  insignificant  Be  indiffusion.  Samples  implanted  to 

higher  Be  concentration  showed  increased  junction  depth.  For  some  reason  which  is  not  yet 

fully  understood,  the  junction  depth  in  samples  implanted  with  a multi-energy  Be  implant  were 

generally  shallower  than  the  junction  depths  in  samples  implanted  with  a single  400-keV  implant. 

From  the  data  it  appears  that  the  maximum  implanted  Be  concentration  for  a flat  profile  without 

18  -3 

significant  indiffusion  is  about  3 x 10  cm  . This  result  is  qualitatively  similar  to  results 

obtained  on  GaAs  (see  Refs.  8 through  11)  and  is  in  fairly  good  agreement  with  recently  obtained 
4 

SIMS  measurements. 

Table  1-1  summarizes  the  sheet  carrier  concentration  and  mobilities  that  have  been  obtained 
on  n-type  InP  and  In,  Ga  As  P,  samples  implanted  with  the  multi-energy  Be  schedule  shown. 

iQ  _ 7 

This  schedule  should  result  in  an  as -implanted  Be  concentration  of  *3  x 10  cm  . For  InP, 
the  highest  sheet  concentration  was  obtained  on  samples  implanted  at  room  temperature  and 
annealed  at  750°C.  As  noted  above,  these  samples  showed  insignificant  indiffusion  compared 
with  samples  implanted  at  150 °C  which  had  a lower  sheet  concentration.  These  results  may 
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TABLE  1-1 


MULTI-ENERGY  BERYLLIUM  IMPLANTS  (FLAT  PROFILE  *3  X 1018  Be  cm"3) 


IN  InP  AND  In.  Go  As  P, 

I -x  x y 1-y 


Material 

Implant 

Temperature 

CO 

Anneal 

Temperature 

CO 

N 

s 

(cm  2) 

»*. 

(cm /''V-iec) 

InP 

25 

700 

1.3X  1014 

62 

InP 

25 

750 

2.2  X 1014 

66 

InP 

150 

750 

1 .8  X 1014 

87 

’0.24As0.52P0.48 

25 

700 

1 .2  X 1014 

82 

j0.25As0.52P0.48 

25 

700 

1 .5  X 1014 

81 

i i o i AO 

All  samples  implanted  with  1.8X10  cm  at  400  keV,  1.2X10  cm  at  200  keV, 
7. OX  1013  cm'2  at  100  keV,  and  4.0X  lo'3  cm'2  at  50  keV. 


indicate  that  outdiffusion  to  the  surface  also  takes  place  when  indiffusion  is  observed,  as  has 
been  shown  to  take  place  in  GaAs  (see  Ref.  10).  Stripping  measurements  on  a sample  implanted 
at  room  temperature  and  annealed  at  750°C  indicate  that  the  hole  concentration  is  fairly  uni- 
form  in  the  flat  portion  of  the  implanted  region,  with  p « (1.5  to  2.0)  x 10  cm  . Although  the 
sheet  concentrations  obtained  in  quaternary  samples  which  are  annealed  at  700  °C  are  not  as 
high  as  the  best  observed  in  InP,  they  are  comparable  to  that  obtained  in  the  InP  samples  an- 
nealed at  700°C.  „ _ _ ,,  _ . _ 


J.  P.  Donnelly  G.  A.  Ferrante 
C.  A.  Armiento  S.  H.  Groves 


B.  P-N  JUNCTION  DIODES  IN  InP  AND  In  Ga  As  P,  FABRICATED 
BY  BERYLLIUM-ION  IMPLANTATION  1-x  x y *'y 


In  this  section  we  report  results  of  p'-n  junction  diodes  formed  by  implantation  of  Be  in 
n-type  InP  and  In^ ^Ga^^As^Pj . Mesa  and  planar  InP  diodes,  which  exhibit  low  leakage  cur- 
rents and  abrupt  breakdowns,  have  been  fabricated.  In  addition.  In,  Ga  As  P,  (E  * 1.0  eV) 

1-x  x y 1-y  g 


mesa  diodes  have  been  made  which  exhibit  uniform  breakdown  over  the  area  of  the  device. 


A / -J 

InP  samples  were  cut  from  (100)  n-type  (n  « 2 x 10  cm'  ) pulled  single  crystals. 


Prior  to  implantation,  these  samples  were  polished  and  etched3  in  a 1 :1 :5:1  mixture  of 
HAc:HC10^:HN03:HC1.  Inj_xGaxAs  Pj  samples  were  n-type  [n  « (0.5  to  1.0)  x 1016  cm'3] 
epitaxial  layers  grown  by  liquid-phase  epitaxy  on  (100)  n+  InP  substrates.6  Beryllium  was  im- 


planted at  room  temperature  into  the  etched  InP  surfaces  and  into  the  as-grown  surfaces  of  the 
In1_xGaxASyP1_y  layers.  Details  of  the  implantation  and  anneal  procedure  can  be  found  in 
Sec.  A above  and  in  Ref.  7. 

Mesa  diodes  were  fabricated  by  implanting  the  entire  surface  of  the  sample  and  subse- 
quently etching  mesas  using  a 1 -percent  bromine -methanol  solution.  Planar  diodes  were 
fabricated  by  implanting  through  holes  opened  in  a 5-p.m-thick  layer  of  photoresist  and  a 
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2000-A  layer  of  pyrolytic  Si02<  In  both  structures,  ohmic  contacts  to  the  p-type  Be-implanted 
layers  were  made  using  microalloyed  evaporated  Au-Mg.  Contacts  to  the  back  of  the  n-type 
substrates  were  made  using  microalloyed  Au-Sn. 

Initial  diodes  were  fabricated  using  the  following  multi-energy  Be  implant  schedule: 

1.8  x 1014  cm*2  at  400  keV,  1.2  x 1014  cm'2  at  200  keV,  7 x 1013  cm'2  at  100  keV,  and 

13  -2  12  13 

4x10  cm  at  50  keV.  According  to  LSS  range  theory,  * this  implant  schedule  should 
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yield  a fairly  uniform  as-implanted  Be  concentration  of  3 x 10  cm  and  a junction  depth  of 
2 pm.  Measurement  of  the  hole  concentration  in  the  implanted  region  and  the  junction  depth 
have  shown  that  this  implant  schedule  results  in  insignificant  indiffusion  of  the  implanted  Be 
(see  Sec.  A above). 

Figures  I-4(a)  and  (b)  show  the  I-V  characteristics  of  a 10-mil-diam  mesa  and  a 
15-mil-diam  planar  InP  diode,  respectively,  fabricated  using  the  above  implant  schedule. 
Leakage  current  in  both  diodes  is  in  the  subnanoampere  range  for  reverse  biases  out  to  more 
than  10  V.  Near  breakdown,  the  etched  mesa  diodes  exhibited  significantly  more  leakage  cur- 
rent than  the  planar  diode.  A possible  explanation  (but  not  necessarily  the  correct  one)  is  that 
the  planar  diodes  had  a pyrolytic  SiO.,  layer  over  the  exposed  junction  which  may  have  provided 
some  surface  passivation. 

Similar  mesa  diodes  fabricated  in  a 6-pm-thick  layer  of  In^  75GaQ  25AsQ  52 P0  48 
(n  a 8 x 1015  cm*3)  exhibited  breakdown  voltage  as  high  as  92  V.  Figure  1-5  shows  the  effec- 
tive concentration  vs  depletion  width  of  such  an  InQ  75GaQ  -^ASq  52P0  48  diode  obtained  from 
capacitance-voltage  measurements.  Using  the  abrupt  depletion  model,  the  effective  concentra- 
tion (Neff)  is  equal  to  /(NA  + N^),  where  and  are  the  net  ionized  acceptor  and 

donor  concentrations,  respectively,  at  the  edge  of  the  depletion  region.  In  the  situation  where 
N.  » N„,  N vs  depletion  width  is  essentially  a plot  of  net  donor  concentration  as  a function 
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of  distance.  Figure  1-5  indicates  that  the  p-n  junction  is  fairly  abrupt.  Although  an  exact  fit 
has  not  been  obtained,  the  slight  grading  that  is  observed  is  about  what  is  expected  from  the 
as-implanted  Be  profile  (which  to  first  order  falls  off  with  a Gaussian  shape). 

To  determine  if  breakdown  in  the  etched-mesa  InP  and  InQ  75GaQ  25ASq  52Pq  48  d*odes 
was  due  to  uniform  avalanche  breakdown,  scanned  photoresponse  measurements  were  performed 
using  the  6328-A  line  of  a He-Ne  laser.  Although  some  of  these  diodes  exhibited  uniform  photo- 
v current  gains  of  up  to  3,  many  diodes  showed  edge  breakdown  and/or  an  enhanced  edge  response. 

As  might  be  expected,  inspection  of  the  mesas  indicated  that  the  slope  of  the  mesa  at  the  junc- 
tion generally  determined  whether  or  not  a diode  showed  these  edge  effects. 

To  eliminate  these  edge  effects  so  that  the  intrinsic  uniformity  of  the  Be-implanted  junctions 
could  be  determined,  diodes  with  shallower  junctions  were  made  in  a thin  In1_xGaxASyP^_y 
layer  which  was  "punched-through"  or  fully  depleted  before  breakdown.  The  shallow  junctions 

permit  more-reproducible  etching  of  mesas  with  near-vertical  walls  at  the  junction,  which 
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minimize  edge  response  and  field  crowding.  Punch- through  operation  further  minimizes  the 

effect  of  any  residual  field  crowding. 

The  In.  ,,Ga.  ,,As.  C,P„  layer  used  for  these  diodes  was  4 pm  thick  and  was 
U.75  0 •cd  U.  dc  U.4o  2 43  _2 

implanted  with  Be  doses  of  6 x 10  cm"  at  70  keV  and  5 x 10  cm  at  50  keV.  Accord- 
ing to  LSS  range  statistics,  this  implant  schedule  should  yield  a junction  depth  of  0.6  pm. 
However,  cleaved,  stained  mesas  indicate  that  some  diffusion  of  Be  has  occurred,  resulting 
in  a junction  depth  of  1.2  to  1.5  pm. 


EFFECTIVE  CONCENTRATION  (cm 


(o)  (b) 

Fig.  1-4.  Current-voltage  characteristics  of  Be-implanted  InP  diodes: 
(a)  10-mil-diam  mesa  diode,  (b)  15-mil-diam  planar  diode. 


Fig.  1-5.  Effective  concentration  as 
a function  of  depletion  width  for  an 

In0.75Ga0.25As0.52P0.48  mesa  diode< 


Fig.  1-6.  Effective  concentration  as 
a function  of  depletion  width  for 

an  In0.75Ga0.25As0. 52P0.48  Punch- 
through  mesa  diode. 
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Figure  1-6  shows  the  effective  concentration  vs  depletion  width  obtained  from  capacitance- 
voltage  measurements  on  a 10-mil-diam  etched-mesa  diode.  Even  though  there  was  apparently 
some  Be  diffusion,  these  measurements  indicate  that  the  junction  is  still  abrupt.  The  sudden 
increase  in  concentration  at  2.5  pm  indicates  that  the  depletion  region  has  reached  through  to 
the  n -InP  substrate.  This  punch-through  occurs  at  45  V,  while  breakdown  occurs  at  72  V. 

Figures  I-7(a)  and  (b)  show  the  scanned  photoresponse  of  one  of  these  devices  at  a reverse 
bias  of  69  V,  which  was  taken  using  the  1.15-pm  line  of  a He-Ne  laser.  Similar  results  were 
obtained  using  the  0.6328-pm  line.  Figure  I- 7 (a ) shows  the  photoresponse  superimposed  on  the 
raster  scan  pattern,  while  (b)  shows  the  horizontal  scans  superimposed  on  each  other  so  that 
fluctuations  in  photoresponse  are  more  easily  detectable. 

These  scans  indicate  a uniform  photoresponse  gain  of  2 (compared  with  the  response  at 
50  V)  over  the  entire  area  of  the  device,  without  any  indication  of  edge  response  or  edge  break- 
down. The  depression  in  the  photoresponse  in  the  center  of  the  device  is  due  to  the  Au-Mg 


(o) 


<k> 


Fig.  t-7.  Photoresponse  scans  at  1.15  pm  of  an  In^  ^Ga^  52^0  48  Punc^~ 

through  mesa  diode:  (a)  3-dimensional  display,  (b)  2-dimensional  display. 
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contact  and  the  lead  wire.  These  measurements  indicate  that  junctions  exhibiting  uniform  ava- 
lanche breakdown  can  be  fabricated  by  Be-ion  implantation. 

C.  A.  Armiento 
J.  P.  Donnelly 
S.  H.  Groves 

C.  GROWTH  OF  NOMINALLY  UNDOPED  InP  AND  InQ  gGa0  2AsQ  gP0  g ALLOYS 

Two  simple  and  reproducible  liquid-phase  epitaxial  growth  techniques  have  been  used  to 
grow  InP  and  InGaAsP  alloys  with  NQ  - NA  at  the  low  1015  cm'3  level.  This  represents  nearly 
an  order-of-magnitude  reduction  in  net  impurity  concentration  below  that  in  growths  made  with- 
out "in  situ"  purification.  The  first  growth  technique  uses  ~5  x 10'7  atm  of  H20  in  the  H2 
growth  atmosphere  to  convert  most  of  the  Si  in  the  growth  solution  to  Si02<  The  second  tech- 
nique uses  PH3  as  a source  of  P,  but  a bake  prior  to  growth  is  necessary.  Results  similar  to 
ours  for  this  technique  were  reported  earlier15’16  but  without  an  explanation  of  the  purification 
mechanism.  We  have  gathered  evidence  which  suggests  that  H20  is  inadvertently  added  with 
the  PII3,  and  the  purification  in  this  process  is  accomplished  by  the  same  mechanism  as  in  the 
H20-addition  process. 

Growth  of  InP  with  77  K electron  mobilities  in  the  40,000  to  60,000  cm2/V-sec  range,  with 
^D^A  between  2.5  and  6,  is  readily  achieved  by  either  process.  For  In^  gGa^  2As  P 
the  77  K mobilities  are  12,000  to  14,000  cm2/V-3ec  with  ND/NA<  2.  The  H20-addition  tech- 
nique offers  greater  control  and  has  given  the  lowest  values  of  ND  - NA>  On  the  other  hand, 
the  PH3 -source  technique  offers  the  attractive  side  benefit  of  providing  a protective  partial 
atmosphere  of  phosphorus  to  prevent  thermal  etching  of  InP  substrates.  A fuller  discussion  of 
this  work  has  been  accepted  for  publication.17 

S.  H.  Groves 
M.  C.  Plonko 

D.  TRANSPORT  PROPERTIES  AND  IMPURITY  LEVELS  OF  InP 

AND  In0.8Gao.2As0.5P0.5 

Hall-coefficient  and  resistivity  measurements  as  a function  of  temperature  have  been  made 
on  InP  and  In0.8Ga0.2As0.5P0.S  aUoys  grown  by  liquid-phase  epitaxy  (LPE).  These  materials 
were  grown  with  one  of  the  "in  situ"  purification  procedures  discussed  in  Sec.  C above,  and  they 
have  net  shallow  impurity  concentrations  Np  - NA  in  the  low  1015  cm'3  range.  Low -magnetic- 
field  measurements  were  made  at  H = 5 kG  with  an  automatic  measurement  apparatus16  for 
temperatures  between  10  and  300  K.  The  effective  electron  concentration  n was  determined 
from  the  measured  Hall  coefficient  R by  assuming  that  the  Hall-coefficient  scattering  factor 
r is  unity:  n = r/Re  = 1/Re,  where  e is  the  electronic  charge.  These  data  are  shown  by  dots 
on  the  main  plots  in  Fig.  1-8  for  InP  and  in  Fig.  1-9  for  an  InQ  gGa0  2AsQ  5PQ  g alloy. 

A fit  of  the  extrinsic  carrier  activation  to  the  data  of  Figs.  I-8  and  1-9  has  been  accom- 
plished by  first  assuming  a single-donor  species,  and  second,  by  ignoring  excited  impurity 
states,  as  is  appropriate  to  the  range  of  purity  considered  here.19  With  these  assumptions,  the 
unknown  parameters  that  determine  the  electron  concentration  are  Nq,  Na,  and  E^  (the  donor 
and  acceptor  concentrations  and  the  donor  activation  energy,  respectively).  The  quantity 
ND  ~ NA  was  estimate<J  from  the  high-temperature  n values  in  Figs.  1-8  and  1-9,  and  with  this 
the  separate  values  of  NQ  and  NA  were  found  from  a mobility  analysis,19  which  is  based  on  the 
assumption  that  at  some  temperature  the  mobility  is  entirely  limited  by  ionized  impurity 
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Fig.  1-8.  Effective  electron  carrier  con- 
centration vs  1000/T  for  InP.  Solid  curve 
is  calculated  using  Ep  = 3.7  meV.  Insert 

shows  high-magnetic -field  data;  curve  is 
calculated  using  a deep-level  concentration 

of  5 x 10*^  cm  ^ and  energy  of  0.29  eV. 


Fig.  1-9.  Effective  electron  carrier  concen- 
tration vs  1000/T  for  Inp  gGa^  2^so  5Po  5‘ 
Solid  curve  is  calculated  using  ED  = 0.  In- 
sert shows  high-magnetic-field  data;  curve 
is  calculated  using  a deep-level  concentra- 
14  -3 

tion  of  3 X 10  cm  and  energy  of  0.12  eV. 


scattering.  The  conduction-band  effective  mass  m * and  static  dielectric  constant  e used  for 

co 

this  analysis  are  m*  = 0.79  m and  e = 12.4  for  InP,  and  m*  = 0.061  m and  e = 13.2  for 
c 00  c 00 

InQ  gGaQ  2Aso  5Po  5*  where  for  the  all°y  is  taken  from  the  magnetoatoorption  results  pre- 
sented in  Sec.  F below,  and  eQ  for  the  alloy  is  estimated  from  an  interpolation  formula  using 
the  reported  values  of  cQ  for  GaAs,  InP,  InAs,  and  GaP.  This  procedure  yielded  the  values 
Nq  = 3.7  x lO*'’  cm  and  = 6.0  x 10*^  cm'^  for  the  InP  sample  from  the  fit  to  the  mobility 
at  T =<  40  K,  and  N^  = 6.0  X 10*  cm  ^ and  = 3.0  X 10*^  cm  ^ for  the  alloy  from  the  fit  to  the 

mobility  at  T » 20  K.  With  these  values  of  and  in  the  charge  neutrality  equation,  E^  was 

adjusted  to  give  an  acceptable  fit  of  calculated  to  measured  electron  concentration.  The  curves 
in  the  main  portion  of  Figs.  1-8  and  1-9  give  the  calculated  curve  for  InP  with  E^  = 0.0037  eV, 
and  for  In^  gGa^  2Aso  5P0  5 ejj  “ respectively.  Boltzmann  statistics  are  adequate  for 
the  InP  case,  but  Fermi-Dirac  statistics  are  necessary  for  the  charge  neutrality  calculation  in 
the  case  of  the  alloy.  The  deviation  of  the  experimental  points  from  the  calculated  curve  at  low 
temperatures,  especially  apparent  in  Fig.  1-8,  is  due  to  the  onset  of  a competing  conduction 
process  such  as  impurity  band  conduction. 

The  behavior  of  the  real  electron  concentration  with  temperature  above  100  K for  both  the 
InP  and  InGaAsP  alloys  is  obscured  by  the  temperature  dependence  of  the  Hall-coefficient 
scattering  factor  due  to  the  increasing  importance  of  polar  mode  phonon  scattering  above  that 
temperature,  an  effect  that  has  been  studied  in  the  case  of  high-purity  GaAs  (see  Ref.  20).  To 
separate  the  carrier  activation  from  scattering  effects  it  is  desirable  to  measure  the  Hall  coeffi- 
cient in  the  high -magnetic -field  limit,  in  which  the  scattering  factor  is  unity.  Measurements 
have  been  taken  at  the  National  Magnet  Laboratory  with  H =*  150  kG,  well  into  the  high-field 
regime,  for  temperatures  between  100  and  400  K,  and  the  resulting  data  are  plotted  in  the 
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inserts  of  Figs.  1-8  and  1-9.  An  effect  that  is  difficult  to  see  on  the  scale  of  these  figures,  but 
is  present  in  all  five  samples  measured,  is  the  extra  activation  of  electrons  discernible  at 
T > 250  K.  To  fit  this,  we  have  added  a deep  level  of  unknown  energy  and  concentration  to  the 
j charge  balance  equation.  After  some  trial  and  error,  we  realized  that  deep  levels  giving  a 

peak  in  the  photoluminescence  (see  Sec.  E below)  for  both  InP  and  InGaAsP  al'oys  were  of  the 
correct  energies  to  produce  the  extra  carrier  activation  if  these  levels  were  taken  as  lying 
closer  to  the  conduction  band  than  the  valence  band.  With  the  energy  difference  between  near- 
band-gap luminescence  and  the  deep-level  luminescence,  0.29  eV  for  InP  and  0.12  eV  lor 
InGaAsP,  taken  as  the  deep-level  ionization  energies,  we  can  generate  the  curves  shown  in  the 

inserts  of  Figs.  1-8  and  1-9  with  deep-level  concentrations  of  5 x 1015  cm'3  for  InP  and 
a 3 

3X10  cm  for  the  InGaAsP  alloy.  The  reduced  number  of  these  levels  in  the  alloy  com- 
pared with  InP  is  in  qualitative  agreement  with  the  intensities  of  the  photoluminescence  lines, 
shown  in  Sec.  E below. 


An  interesting  possibility  is  that  the  deep  level  is  due  to  oxygen.  Several  workers  growing 
epitaxial  InP  have  noticed  that  growth  in  an  atmosphere  with  more  Oz  or  HzO  than  usual  causes 
an  increase  in  the  ratio  of  electron  concentration  at  300  K to  that  at  77  K,  n3o0/n77-  This 
effect  may  be  caused  by  an  increased  number  of  these  deep  levels,  although  it  is  recognized 
that  the  concentration  and  ionization  energy  of  the  shallow  donors  as  well  as  the  compensation  - 
factors  that  may  vary  with  the  02  or  H20  concentration  in  the  growth  tube  atmosphere  - also 
affect  the  quantity  n /n,_. 

300'  77  S.  H.  Groves 


E.  PHOTOLUMINESCENCE  OF  InP  AND  In0>gGa0  2AsQ  5PQ  5 ALLOYS 

A preliminary  study  has  been  made  of  photoluminescence  from  nominally  undoped  InP  and 
In0.8Ga0.2As0.5P0.5  alloYs  grown  by  liquid-phase  epitaxy  (LPE).  This  study  has  been  done 
with  the  samples  at  77  K and  with  an  apparatus  selected  for  ease  of  operation  over  a broad 
range  of  photon  energies,  0.3  to  1.6  eV,  at  some  sacrifice  of  resolution.  A krypton  laser  has 
provided  photoexcitation  of  0.5  W at  0.6741  pm  wavelength.  The  luminescent  radiation  has  been 
passed  through  a prism  monochromator  and  imaged  on  a thermoelectrically  cooled  PbS  detector. 

Figure  1-10  shows  the  luminescence  intensity  vs  photon  energy  for  a sample  with  electrical 
characteristics  ND  - NA  = 2.8  x 1015  cm'3  and  = 3.7  x 104  cm2/V-sec,  where  ND  - NA  is 
the  net  impurity  concentration  and  p??  is  the  electron  mobility  at  77  K.  The  peak  at  ~1.4  eV 
involves  transitions  to  shallow  donors  and  acceptors  and,  with  higher  resolution,  would  sepa- 
rate into:  (1)  a sharp,  high-intensity  peak  at  1.4  eV  due  to  band-to-band  and  shallow-donor-to- 
valence-band  transitions,  (2)  a weaker  peak  35  to  40  meV  below  the  1.4-eV  peak  and  due  to 
conduction-band-to-shallow -acceptor  transitions,  and  (3)  possibly  a transition  75  to  80  meV 
below  the  1.4-eV  peak  and  very  much  weaker.  These  details  can  be  seen  in  the  spectra  re- 
ported by  workers  at  the  Royal  Radar  Establishment  ^ for  LPE  grown  samples  measured  at 
77  K and  in  work  we  have  done22  in  an  apparatus  with  higher  resolution  and  with  the  samples 
near  4 K.  Very-high-resolution  near-gap  luminescence  has  been  reported  for  LPE  InP  samples 
at  1.8  K by  Hess  et  al.23 

The  second  peak  in  Fig.  1-10,  at  1.12  eV,  is  probably  of  the  same  origin  as  the  broad  peak 
at  1.17  eV  in  Ref.  21.  This  was  seen  in  most,  but  not  all,  of  the  undoped  LPE  material  investi- 
gated in  that  study  and  was  seen  in  the  bottom  portion  of  their  undoped,  melt -grown  InP  crystals. 
The  hypothesis  was  made  in  Ref.  21  that  this  line  comes  from  impurity/P-vacancy  complexes. 
From  our  analysis  of  extrinsic  carrier  activation  discussed  in  Sec.  D above,  we  believe  that  this 
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transition  involves  a level  located  ~0.3  eV  below  the  conduction  band.  As  mentioned  above,  an 

investigation  is  under  way  to  test  the  hypothesis  that  this  level  is  due  to  oxygen. 

The  lowest  energy  peak  in  Fig.  1-10  occurs  at  0.49  eV.  Later  runs  with  this  apparatus  were 

extended  to  0.3  eV,  and  all  samples  showed,  in  addition,  a similar  type  of  transition  at  0.35  eV. 

24 

On  the  basis  of  an  electron  spin  resonance  analysis,  the  0.35-eV  line  has  been  assigned  to  am 
intra-center  transition  between  the  five  and  five  E levels  of  the  Fe2  + (3d^)  impurity.  The 
origin  of  the  0.49 -eV  level  is  not  known  but  may  be  radiation  from  a higher  excited  level  of  the 
Fe2  + manifold.  Our  undoped  layers,  5 to  10  pm  thick,  were  grown  on  the  Fe-doped  InP  sub- 
strates. However,  the  profile  of  the  Fe  in  the  sample  and  the  source  of  the  Fe  giving  rise  to 
the  photoluminescence  lines  have  not  been  determined. 
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Fig.  I- 1 0.  Photoluminescence  spectrum  Fig.  I— 11.  Photoluminescence 

for  undoped  InP  grown  by  LPE.  spectrum  for  an  undoped  alloy 

of  In0.8Ga0.2As0.5P0.5  grOWn 
by  LPE. 

Figure  1-11  shows  the  luminescence  spectrum  for  an  InQ  8GaQ  -,AsQ  5PQ  5 alloy  with  elec- 
trical characteristics  = 5.8  x 1015  cm  3 and  p77  = 1.2  x 10^  cm2/V-sec.  The  relative 

intensity  scales  in  Figs.  1-10  and  1-11  are  not  necessarily  the  same,  but  probably  are  within  a 
factor  of  2 of  each  other.  The  main  difference  between  the  traces  in  Figs.  1-10  and  I- 1 1 , other 
than  the  energy  shift  due  to  the  smaller  bandgap  of  the  alloy,  is  the  reduced  intensity  in  the  sec- 
ond peak  of  Fig.  I— 11.  Again,  from  the  extrinsic  activation  analyses  above,  we  attribute  this 
line  to  a level  located  0.12  eV  below  the  conduction  band.  The  reduced  concentration  of  this  level 
in  the  alloy  compared  with  that  in  InP  is  in  qualitative  agreement  with  the  results  of  Sec.  D above. 

S.  H.  Groves 
J.  G.  Mavroides 
D.  F.  Kolesar 
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Fig.  1-12.  Voigt  configuration  magneto-transmission  spectrum.  Bars  at  bottom 
indicate  theoretical  energies,  and  their  heights  are  proportional  to  square  of 
matrix  elements. 


i u h ii  [100]  /ry  /y 

• EXPERIMENT  ///  yy  ~ 

THEORY  ( heavy  holo)  /jf/  /y  /; 

THEORY  (light  hole*)  ///  Z' 


Fig.  1-13.  Transmission  minima  vs 
magnetic  field.  Theoretical  curves 
include  estimated  excitation  binding 
energy,  which  at  H = 0 is  3 meV. 
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F.  INTERBAND  MAGNETOOPTICAL  MEASUREMENTS 

ON  Rij  _xGaxAsyPly  ALLOYS 

Oscillatory  interband  magneto-transmission  in  the  1.2-nm  wavelength  region  has  been  mea- 
sured on  samples  of  In^  7yGaQ  23^ao  52^*0  48  8rown  by  liquid-phase  epitaxy  on  InP  substrates. 
Figure  1-12  shows  the  oscillatory  transmission  made  in  the  Voigt  configuration,  E ||  H ||  [100], 
for  H = 130.5  kG,  where  E is  the  optical  electrical  field  vector.  The  arrows  show  the  energies 
and  strengths  of  transitions  predicted  by  the  theory  with  band  parameters  determined  from  this 
study.  Figure  1-13  shows  the  spectrum  generated  by  plotting  energies  of  absorption  minima 
vs  H.  The  band  parameters  determined  by  the  fit  of  theory  to  the  experimental  data  are 
Eg  = 1.07  eV,  m*  = 0.061  mQ,  and  Ep  = 17.6  eV,  where  these  parameters  are,  respectively, 
the  direct  energy  gap  at  T » 20  K,  the  conduction-band  effective  mass,  and  the  It  • fT  interaction 
energy,  proportional  to  the  square  of  the  momentum  matrix  element.  We  are  presently  extend- 
ing the  magnetooptical  measurements  both  to  other  InGaAsP  alloys  lattice -matched  to  InP  and  to 
such  techniques  as  electric -field -modulated  magnetoreflection.  A more-detailed  report  of  this 


work  has  been  accepted  for  publication. 


K.  Alavit 

R.  L.  Aggarwalt 

S.  H.  Groves 


t Francis  Bitter  National  Magnet  Laboratory,  M.I.T. 
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II.  QUANTUM  ELECTRONICS 


A.  EMISSION  CROSS  SECTION  AND  FLASHLAMP-EXCITED 

NdP5014  LASER  AT  1.32  pm 

1-3 

Stoichiometric  neodymium  materials  have  been  the  subject  of  considerable  investigation 
based  on  their  potential  application  in  efficient  miniature  solid  state  lasers.  Their  high  Nd  con- 
centration causes  efficient  pump  absorption,  yielding  high  gain  per  unit  volume,  and  their  free- 
dom from  significant  concentration  quenching  results  in  long  fluorescent  lifetimes.  We  have 
measured  the  emission  cross  section  and  the  room-temperature  flashlamp-pumped  laser  action 

of  NdP-O . . at  1.32  pm.  Previous  laser  experiments  on  similar  materials  are  either  at  the 

^ 4 + + 

conventional  1.05-pm  line  or  have  used  Ar  (Ref.  5)  or  Kr  (Ref.  6)  lasers  as  excitation  sources. 

We  also  compared  the  1.32-pm  laser  performance  with  that  of  the  conventional  1.05-pm  line. 

The  effective  emission  cross  section  determined  from  fluorescence  measurements  was  in  good 
agreement  with  the  value  found  from  laser-threshold  data.  The  single-crystal  1-mm2  x 7-mm 
NdP5014  laser  rod  delivers  a few  millijoules  multimode  output  energy  per  pulse  at  1.32  pm, 
and  the  threshold  is  as  low  as  180  mj  of  electrical  input  energy  with  a 0.4-percent  transmitting 
output  mirror.  The  1.32 -pm  wavelength  is  attractive  for  optical  fiber  communications  since  at 
this  wavelength  transmission  loss  of  silica  fibers  is  low  and  the  material  dispersion  is  favor- 
able.7 The  lasing  resonant  losses  at  1.32  pm  are  negligible  compared  with  those  at  1.05  pm, 
since  the  4I13/2  lower  laser  level  of  the  1.32-pm  line  is  about  3900  cm”1  above  ground  level. 
Other  advantages  compared  with  1.05-pm  operation  are  larger  mode  volumes  at  1.32  pm  and 
lower  Rayleigh  scattering  loss.  Compared  with  quaternary  semiconductor  diode  lasers  at  this 
wavelength,  stoichiometric  Nd  lasers  might  be  a feasible  alternative  because  of  smaller  beam 
divergence,  narrower  laser  linewidth,  and  higher  output-pulse  energies.  Relative  to  its 
Nd:  YAG  counterpart,  the  wider  linewidth  of  NdP5C>14  would  have  shorter  modelocked  pulse 
capabilities. 

Single  crystals  of  NdPg014  were  grown  from  hot  phosphoric  acids  in  which  the  rare-earth 
oxide  had  been  dissolved.  Details  of  the  crystal  growth  were  reported  previously.9 

Fluorescence  measurements  were  performed  on  NdPg0^4  samples  excited  by  a chopped 
argon-ion  laser  at  5145  A.  The  excitation  radiation  was  focused  on  a 1.5-mm  "c"  plate  of 

NdP-O.  . of  1-mm2  cross-section  area  placed  close  to  the  entrance  slit  of  the  0.25-m  Jarrell- 

5 14 

Ash  spectrometer  with  a 590  grooves/mm  grating  blazed  at  1.2  fim.  After  transmission  through 
a long-pass  filter,  infrared  fluorescence  was  detected  by  a room-temperature  Ge  photodiode  of 
400-MHz  bandwidth  and  a lock-in  amplifier.  An  unpolarized  fluorescence  spectrum  at  1.32  pm 
is  shown  in  Fig.  II—  1. 

Overall  system  calibration  of  detector  and  spectrometer  was  carried  out  with  a tungsten 
lamp,  using  the  appropriate  emissivity  and  blackbody  curves.  The  calibration  was  further 
checked  against  the  known  specifications  of  the  detector  spectral  response  and  the  unpolarized 
grating  efficiency.  The  effective  emission  cross  section  at  1.32-pm  transition  was  determined 
by  comparing  the  peak  emission  intensity  of  the  1.32-pm  line  with  that  at  1.05  pm,  whose  cross 
section  is  already  known.  The  relation  between  the  cross  sections  of  the  two  transitions  is 
approximately1  ° 


As  fluorescence  data  are  compared  under  identical  excitation  conditions,  reabsorption  correc- 

1 1 

tion,  which  is  of  importance  for  thick  samples,  was  canceled  out  in  the  relative  measurement. 
With  the  system  calibration  correction,  a ratio  of  the  1.05-  to  1.32-pm  peak  cross  sections 
u./a2  = 5.6  was  obtained  from  the  fluorescence  data.  By  using  a value  of®,  n,  = 2.5  x 10  cm 
at  1.05  pm  (Ref.  12)  the  cross  section  at  1.32  pm  was  found  to  be  33  = 4-5  * 10  cm  . 

Laser  rods  of  NdP-O. . were  fabricated  with  "c"  faces  polished  flat  (X/10)  and  parallel  to 

3 3+  4 

10  sec,  and  "a"  and  "b"  faces  fine  ground.  The  fluorescence  lifetime  of  the  Nd  F3/2  mani_ 

fold  was  135  psec  as  measured  by  a laser  pumped  dye  laser.  The  uncoated  laser  rod  was 

placed  inside  an  aluminum-coated  reflector  close  to  a Xe  flashlamp  which  was  driven  by  a single 

4 

mesh  circuit  of  10-pF  capacitance  and  7 5-pH  inductance  described  previously.  (The  mode  waists 
of  the  10-cm  confocal  cavity  are  129  and  183  pm  at  the  crystal  and  mirrors,  respectively.)  The 
laser  mirrors  (Optasil)  should  have  negligible  absorption  at  1.32  pm.  Output  mirrors  of  99.6-, 

99. 0-,  and  97.9-percent  reflectivities  at  1.32  pm  were  used  in  conjunction  with  a highly  reflect- 
ing back  mirror.  All  the  ditlectric  coatings  had  >85-percent  transmission  at  1.05  pm.  With 
the  99.0-percent  reflecting  mirror,  an  output  energy  of  1.5  mj  per  pulse  at  a 1-Hz  repetition 
rate  was  measured  with  a power  meter.  The  laser  linewidth  at  1.32  pm  was  14  cm-1  at  2 -percent 
output  coupling  as  measured  by  the  Ge  photodiode  at  the  output  slit  of  a 0.25-m  Jarrell-Ash  spec- 
trometer. As  in  the  case  of  1.05-pm  operation,  the  1.32-pm  laser  polarization  for  the  "c"-axis 
rod  was  along  the  "b"  direction.  With  the  1-mm2  x 7 -mm  laser  rod,  laser  performance  at 
1.32  pm  was  compared  with  1.05-pm  operation  under  identical  cavity  configurations.  With  an 
output  coupling  of  ~2  percent,  the  threshold  energy  at  1.32  pm  was  340  mj;  with  2.5-percent 
coupling,  that  at  1.05  pm  was  160  mj.  At  this  coupling  level  the  multimode  output  energy  per 
pulse  was  2.7  and  8.2  mj,  respectively,  with  1.6  J of  electrical  energy  input.  The  slope  effi- 
ciencies were  0.24  and  0.52  percent  at  1.32  and  1.05  pm,  as  shown  in  Figs.  II— 2 and  II— 3-  Total 


internal  cavity  losses  were  estimated  at  2 percent  in  both  cases,  with  the  1.05-pm  operation 
having  an  additional  3-percent  resonant  absorption  loss,  while  that  at  1.32  pm  is  negligible. 

From  a comparison  of  the  change  in  threshold  energy  with  varying  mirror  transmission  for  both 
1.05-  and  1.32 -pm  operation,  the  ratio  of  o^/cr.,  = 5*75  was  found.  In  this  case,  the  ratio  refers 
to  the  population-weighted  "b"  polarization  cross  sections  rather  than  the  presumably  mixed 
polarization  ratio  given  above.  Using  only  the  "b"  polarized  cross-section  data,12  a.  nc(b)  = 

— “IQ  2 o 1#U!> 

1.6  X 10  y cm  , giving  a ^ 32(b)  = 2,8  x 10  cm. 

S.  R.  Chinn 
M.  M.  Choy* 

W.  K.  Z wicker* 

B.  Ni:MgF2  LASER  DEVELOPMENT 

We  have  obtained  100  mW  of  CW  output  from  a 3-mirror  cavity  Ni-.MgF,  laser  described 
13  ^ 

earlier,  and  have  observed  extremely  low  crystal  losses.  Repetitive  Q-switching  of  the  laser 

3 

has  yielded  power  enhancements  of  >10  , with  peak  power  outputs  of  20  W. 

The  laser  crystal  used  in  these  experiments  was  doped  with  l-wt%  Ni  and  absorbed  ~7  5 per- 
cent of  the  1.33-pm  pump  radiation.  With  a 1.1 -percent-transmission  output  mirror,  we  found 
a CW  power  output  of  100  mW  for  an  absorbed  pump  power  of  1 W.  Threshold  absorbed  powers 
PTH’  sloPe  efficiencies  A,  and  threshold  wavelengths  *TH<  for  three  different  output-mirror 
transmissions  T,  are  given  in  Table  n-1.  We  have  compared  the  thresholds  and  slope  efficien- 
cies for  the  two  lowest-transmission  mirrors  (for  which  the  same  XTH  was  observed)  and  we 


TABLE  ll-l 

Ni:MgF2  LASER  PERFORMANCE 

VS  OUTPUT  MIRROR  TRANSMISSION 

T 

PTH 

A 

XTH 

(percent) 

(mW) 

(percent) 

(h™) 

0.2 

20 

5.5 

1.67 

1.1 

68 

11.5 

1.67 

6.6 

420 

15.0 

1.62 

have  performed  ringing  frequency  measurements,  all  in  order  to  determine  the  total  loss  in  the 
laser  cavity  excluding  that  due  to  output-mirror  transmission.  Our  best  estimate  for  this  loss 
is  0.3  percent.  The  total  path  length  in  the  laser  crystal  for  one  round  trip  of  the  cavity  is 

2.5  cm.  Thus,  if  all  the  loss  is  distributed  uniformly  in  the  crystal  in  the  form  of  absorption 

-3-1 

and  scattering  losses,  we  find  a loss  coefficient  of  1.2  x 10  cm  . However,  we  feel  that  the 
Ni:MgF2  distributed  loss  is  lower  (perhaps  by  50  percent)  than  this  value  since  the  0.3-percent 
loss  includes  scattering  and  transmission  loss  from  two  mirror  surfaces,  scattering  and  reflec- 
tion loss  from  two  Brewster-angle  surfaces  on  the  Ni:MgF2  laser  crystal  and  two  Brewster  sur- 
faces on  the  fused-silica  Brewster  window,  as  well  as  absorption,  scattering,  and  depolarization 


t Philips  Laboratories,  Briarcliff  Manor,  New  York. 


losses  in  a 1-cm  path  of  fused  silica  in  the  latter  window.  The  ultimate  limit  for  crystal  loss, 
arising  from  phonon-assisted  absorption  processes,  is  roughly  4 x 10"5  cm"1  at  1.67  pm  at  77  K. 

Q-switching  of  the  3-mirror  cavity  Ni:MgF2  laser  was  achieved  using  a 2.5-cm-long  crystal 
of  LiNbOj  cut  for  Brewster-angle  operation  with  optical  propagation  along  the  "z"-axis.  The 
insertion  loss  of  the  Q-switch  was  unexpectedly  high,  about  3 percent,  and  was  presumably  due 
to  strain-induced  birefringence  in  the  LiNbOy  Thus,  with  a 1-percent-transmission  output 
mirror,  only  14  mW  of  CW  output  was  observed  at  1.62  pm  for  an  absorbed  pump  power  of  1 W. 
Roughly  200  V applied  across  the  5-mm-thick  Q-switch  was  sufficient  to  turn  the  laser  off. 
Q-switched  output  pulses  observed  by  a high-speed  HgCdTe  photodiode,  for  a pulse  repetition 
rate  (PRR)  of  100  Hz,  are  shown  in  Fig.  II-4.  Also  shown  in  the  figure  is  the  voltage  applied 
across  the  Q-switch.  The  peak  and  average  power  outputs  at  a 100-Hz  PRR  for  1 W of  absorbed 
pump  power  were  20  W and  2.2  mW,  respectively.  At  1 kHz,  peak  and  average  power  outputs 
changed  to  4 W and  9.5  mW. 


Fig.  H-4.  Q-switched  pulses  from 
Ni:MgF2  laser  along  with  voltage 
(~200  V)  drive  to  LiNbOj  Q-switch. 
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The  power  enhancement  obtained  by  Q-switching,  defined  as  the  ratio  of  Q-switched  peak 

•a 

power  to  CW  power  observed  with  no  Q-switching,  was  1.4  X 10  at  a 100-Hz  PRR.  Computer 
calculations  of  CW-pumped,  repetitively  Q-switched  lasers  predict  for  the  system  discussed 
here  a power  enhancement  of  2.5  x 104.  Some  if  not  all  of  this  discrepancy  is  related  to  the 
transient  elastooptic  effect,  which  is  particularly  noticeable  in  LiNbO,  due  to  the  large  difference 

jr  U -5 

between  static  and  dynamic  electrooptic  coefficients.  ’ Essentially,  the  retardation  of  the 
Q-switch  does  not  drop  to  zero  as  the  voltage  drops  to  zero,  but  oscillates  at  the  mechanical 
resonance  frequencies  of  the  bulk  crystal.  Oscillations  in  crystals  similar  to  the  one  we  used 
have  been  observed  with  periods  of  about  1 psec,  lasting  for  >10  p.sec  (Ref.  16).  Since  our  laser 
system  performance  was  sensitive  to  small  losses,  the  oscillations  in  the  retardation  of  the 
Q-switch  would  noticeably  affect  the  power  enhancement.  An  acoustooptic  Q-switch  using  fused 
silica  as  the  interacting  medium  would  likely  provide  better  performance  than  we  achieved  with 
LiNbOj,  since  the  insertion  loss  could  be  much  smaller,  and  there  would  be  no  problem  with  the 
transient  elastooptic  effect. 
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To  summarize,  we  have  observed  >10 3 power  enhancement  in  the  NiiMgF^  laser  by  repeti- 
tive Q-switching,  though  even  greater  enhancements  may  be  possible  with  a more-optimum 
Q-switch.  The  ability  to  achieve  such  large  enhancements  is  an  important  property  of  the 
Ni:MgF2  laser  for  applications  involving  nonlinear  optical  effects,  such  as  frequency  mixing  and 
Raman  shifting. 

P.  F.  Moulton 
A.  Mooradian 

C.  HIGHER-ORDER  NONLINEAR  PROCESSES  IN  CdGeAs2 

We  report  observation  of  phasematched  fourth-harmonic  generation  (FHG)  in  CdGeAs2  using 
a 0O2  pump  laser  with  a 140 -nsec  pulse  duration.  From  the  measurements,  we  determine  an 
effective  fourth-order  susceptibility  tensor  which  is  much  higher  than  previous  theoretical  esti- 
mates of  the  fourth-order  tensor.47,4** 

Measurements  of  higher-order  nonlinearities  are  complicated  because  of  interference  be- 
tween the  direct  processes  and  cascade  processes  connected  with  lower-order  nonlinearities. 

The  susceptibility  tensors  obtained  by  these  measurements  are  effective  values  that  include  the 
contributions  of  the  cascade  processes.  For  CdGeAs.,,  the  estimated  contributions  of  these 
processes  to  the  effective  third-  and  fourth-order  susceptibility  tensors  are  comparable  to  the 
measured  values. 

Quantitative  studies  of  higher-order  nonlinear  optical  processes  are  important  for  the  under- 
standing of  the  nonlinear  properties  of  dielectric  media.  Ninth-harmonic  generation  has  re- 
cently been  reported  in  a gaseous  sodium-argon  mixture,  and  phasematched  fourth-  and  fifth- 
harmonic  generation  in  solids  has  been  observed  in  lithium  formate20  and  calcite,21  respectively. 
These  experiments  were  conducted  in  the  visible  and  ultraviolet  spectral  regions  with  picosecond 
pump  pulses  from  a modelocked  Nd3+  glass  laser.  The  present  experiment  is  the  first  measure- 
ment of  an  effective  fourth-order  susceptibility  tensor  for  an  infrared  material. 

The  interest  in  CdGeAs2  during  recent  years  derives  from  its  large  second-order  nonlinear 

coefficient.  The  material  is  useful  for  frequency- doubling  of  C02-laser  radiation,  and  energy 

conversion  efficiencies  as  high  as  27  percent  have  been  reported.22  Both  its  second-order23'24 

and  effective  third-order  nonlinear  coefficients  have  been  measured.  Because  it  has  the  chal- 

copyrite  structure  (42 -m  symmetry),  CdGeAs2  has  only  two  independent  fourth-order  nonlinear 

coefficients,  v ' and  \ = v (4)  mef  26\ 

*xyzzz  xxyyyz  Xxxxyz  'n  1 

CdGeAs2  is  one  of  the  few  infrared  nonlinear  materials  with  sufficient  birefringence  to  allow 
phasematched  FHG.  Type  I phasematching  is  satisfied  at  pump  wavelengths  longer  than  9.7  pm, 
and  type  II  phasematching  is  allowed  at  wavelengths  between  12.4  and  15.2  pm  (see  Refs.  17  and 
24).  The  type  I phasematching  condition  can  be  stated  as 

n4°  = nlle)  (II-l) 

where  o and  e refer  to  ordinary  and  extraordinary  waves,  respectively,  and  e is  the  angle  be- 
tween the  wavevector  and  the  optic  axis  (z-axis).  If  0 is  the  azimuth  angle  between  the  x-axis 
and  the  projection  of  the  wavevector  in  the  x-y  plane,27  the  fourth-order  nonlinear  coefficient 
for  type  I phasematching  can  be  written  as 

Xj4*(w.u>,u>,u))  = 2 sin 20  cos  20  [x^zz  sin2e  + 

This  expression  is  maximized  for  an  input  wave  polarized  either 
or  ir/2,  respectively). 


(4)  2 

*xyyy*  COS  0J  ' <n-2> 

in  the  x-z  or  y-z  plane  (0=0 


(4) 

A direct  measurement  of  Xj  (w,u>,u>,a>)  is  not  possible  because  of  interference  between 

the  direct  FHG  process  and  the  two-  and  three-step  cascade  processes  involving  the  lower-order 
(2)  (3) 

nonlinearities  x and  x . The  cascade  processes  are  automatically  phasematched  in  the  same 

direction  [defined  by  Eq.  (II- 1) ] as  the  direct  process  even  though  the  individual  cascade  steps 

28 

are  generally  not  phasematched  in  the  same  direction.  The  effective  nonlinearity  which  is 

1 8 

measured  can  therefore  be  written  as 

x/gjf(4w)  = x/4^(w . w , w , u) ) +b1x(2^(w,w)  x*3^(2u>,u>,w) 

+ b2x*3\w,  w,  u>)  x^Ciu.u)  + b3x^(“.  w)  x^(2“.“>)  X^(3w,u>) 

+ b4  (x(2)(w.u))]2  x(2)(2w.2«)  . (n-3) 

The  last  term  vanishes  in  the  type  I phasematching  geometry  that  maximizes  Xj  (u>,  u),  w,  u>). 

In  this  case  (<p  = x/2),  Eq.  (11-3)  takes  the  specific  form 


*I.4eff(4u)  = -2  8in2e  sin2e  + X^iyz  cos2e  + 


4n2°  [»>)  - n2°]  e<> 


(3) 

sin2  9 + x 

xxzz 

*xxyy 

, (3) 

sin40  + x <3* 

zzzz 

Axxxx 

[X(2)]3 

*xyz 

2n®  (0)  [n®  (9)  - n^(0)]  eQ 


cos40  + 6x^z  sin2e  cos 2 9] 


2n3e(9)  n®  (n^e)  - n.°]  [n*(e)  - n/(e)[  e 2 | 

This  equation  is  valid  as  long  as  the  coherence  lengths  associated  with  the  intermediate  steps 
of  the  cascade  processes  are  shorter  than  both  the  crystal  length  and  the  absorption  length. 

Two  CdGeAs2  crystals  2.7  and  3.8  mm  long,  referred  to  as  A and  B,  were  used  in  the  ex- 
periments. Their  measured  absorption  constants  at  liquid  nitrogen  temperature  at  10.25  and 
2.56  pm  were,  respectively,  a^  = 0.1  cm  4 and  a4  = 4.4  cm-4  for  crystal  A,  and  a^  = 1.9  cm-4 
and  a4  = 9.5  cm  4 for  crystal  B.  Both  crystals  were  cut  for  type  I phasematching.  The  mea- 
sured polarization  of  the  fourth-harmonic  output  beam  was  perpendicular  to  the  input-beam  po- 
larization, as  expected  for  type  I phasematching.  Figure  II- 5 shows  the  measured  wavelength 
dependence  of  the  type  I phasematching  angle  for  FHG.  The  measured  angles  have  an  absolute 

uncertainty  of  ±1".  The  figure  also  shows  a phasematching  curve  calculated  from  refractive 
24 

index  data  taken  at  room  temperature.  Previously,  we  have  shown  that  the  phasematching 

angle  for  second-harmonic  generation  (SHG)  of  CO, -laser  radiation  is  only  weakly  temperature- 

‘ 22 

dependent,  varying  by  <1°  between  room  temperature  and  liquid  nitrogen  temperature.  A 

stronger  temperature  dependence  is  expected  for  FHG  since  the  fourth-harmonic  wavelength  is 

close  to  the  wavelength  corresponding  to  the  CdGeAs,  bandgap,  which  varies  from  2.2  pm  at 

c 29 

room  temperature  to  1.9  pm  at  liquid  nitrogen  temperature.  The  increased  bandgap  at  liquid 
nitrogen  temperature  should  reduce  the  indices  of  refraction  at  the  fourth-harmonic  wavelength. 


PHASEMATCHING  ANGLE  (deg) 


MEASURED  80K 


4 10.« 

GTN  (um) 


TABUE  Cl— 2 

MEASURED  UNEAR  AND  NONLINEAR 
SUSCEPTIBILITIES  OF  CdGeAs2 

X = 12 

Azz'  0 

Refs.  23  and  24 

Xx^/‘o  = 6-1X10'1°m/V 

Refs.  23  and  31 

XlSff/to=7-0)<,0'18m2/v2 

Reft.  17  and  25 

X|(/V«o=7-9X,0'27m3^3 

This  work 

and  consequently  lead  to  a smaller  phasematching  angle.  In  Fig.  II-5,  this  effect  is  particularly 
evident  at  the  shortest  pump  wavelengths  where  the  measured  phasematching  angle  is  ~5°  smaller 
than  calculated.  At  10.6  pm  the  difference  is  3.5°,  which  is  slightly  larger  than  the  difference  of 
2°  for  the  phasematching  angle  of  SHG  with  10.6-pm  radiation. 

The  measured  fourth-harmonic  output  power  at  2.56  pm  vs  input  power  at  10.24  pm  for  the 
two  crystals  with  the  gain  cell  on  and  off  is  shown  in  Fig.  II-6.  Turning  the  gain  cell  off  resulted 
in  partial  modelocking  of  the  laser  pulse,  increasing  the  conversion  efficiency  by  15  to  20  times. 
For  the  modelocked  case.  Fig.  II-6  gives  the  average  peak  power.  The  conversion  efficiency 
saturates  at  power  levels  approaching  the  laser-induced  surface  damage  threshold.  The  satura- 
tion is  less  evident  when  the  gain  cell  is  on.  This  suggests  an  intensity-dependent  effect.  From 
the  measured  input-output  curves  in  the  case  of  no  modelocking,  we  determine  the  effective 
fourth-order  susceptibility  tensor  X^4J{{,  in  mks  units,  as  6.2  x 10-38  and  8.0  x 10-38  m2As/V4 
for  crystals  A and  B,  respectively.  In  order  to  reduce  the  possibility  of  systematic  errors, 

X[ 4ff  for  crystal  A was  also  measured  in  another  setup  using  a different  laser  and  detector. 

The  results  of  the  two  measurements  agreed  within  20  percent.  From  the  experiments  on  both 
crystals,  we  obtain 

2 

(4)  n V <n~38  m As  , .n 

xI,eff  = 7 * 10  — 4 * 40  percent 

for  CdGeAs2. 

The  intrinsic  breakdown  thresholds  in  solids  are  determined  by  electron  avalanche  ioniza- 
30  8 

tion,  which  occurs  at  electric-field  strengths  of  typically  10  V/m.  To  reach  this  intrinsic 

threshold  in  CdGeAs2  probably  requires  C02-laser  pulses  shorter  than  1 nsec.  In  Table  II-2, 
the  measured  Xj  gff/e0  is  compared  with  the  lower-order  nonlinear  susceptibilities.  The  higher- 
order  terms  in  the  expansion  of  the  induced  polarization  do  not  exceed  the  lower-order  terms 
even  at  the  threshold  for  dielectric  breakdown. 

The  interference  terms  for  the  fourth-order  process  can  only  be  estimated,  since  the  sign 
(3) 

and  magnitude  of  each  xiS”  component  are  not  known.  An  upper  bound  can  be  determined  using 
(3) 

the  reported  xeff  value  for  type  II  phasematching,  keeping  in  mind  that  this  value  most  likely  is 
(3) 

greater  than  the  X|yJ  components  because  of  the  interference  effects.  The  estimated  contribution 

to  x/gff/£0  in  Eq.  (II-4)  due  to  the  cascade  terms  is  1.7  x 10-26  m3/V3,  assuming  x^yj « X^/3  K 

xeff  > °*  This  is  the  same  order  of  magnitude  as  the  measured  value  of  Xj  4ffA0* 

The  Xj^(u),  w,  ui , a))A0  for  CdGeAs2  has  been  theoretically  estimated  as  1.3  x 10-2^  m3/V3 
17  -29  3 / 3 

using  a bond-orbital  model,  and  as  0.7  x 10  m /V  from  a simple  anharmonic-oscillator 

model.  These  estimates  are  smaller  than  the  measured  Xj4« A value  by  factors  of  almost 

10  . This  suggests  that  the  interference  terms  make  the  main  contribution  to  Xj  The  theory 

may  be  in  error,  however,  since  theory  and  experiment  differ  by  a factor  of  50  for  lithium  for- 
18 

mate  (even  after  accounting  for  the  inference  effects),  the  only  other  material  for  which 
fourth-order  nonlinear  processes  have  been  studied. 
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D.  HIGH-RESOLUTION  TUNABLE  SUBMILLIMETER  SPECTROSCOPY 

32 

In  a previous  Solid  State  Research  Report,  the  generation  and  detection  of  tunable  sub- 
millimeter  radiation  using  corner-reflector  diode  mixers  was  reported.  Two  source  systems 
were  described:  one  which  generated  tunable  sidebands  by  mixing  a microwave  signal  with  the 
output  of  a submillimeter  laser,  and  one  which  generated  a high-order  harmonic  directly  from 
a microwave  signal.  The  first  system  has  been  developed  into  a tunable  sideband  spectrometer 
for  linear  absorption  spectroscopy  which  has  higher  sensitivity  than  the  harmonic  spectrometer. 

Radiation  from  an  optically  pumped  submillimeter  laser  is  coupled  into  both  the  sideband 
generator  and  the  heterodyne  receiver.  In  the  sideband  generator,  the  laser  radiation  is  mixed 
with  a microwave  signal  of  up  to  45  GHz  in  a corner-reflector-mounted  Schottky-barrier  diode. 
The  antenna  structure  which  couples  the  laser  radiation  into  the  diode  reciprocally  reradiates 
the  sidebands.  A diplexer  separates  the  reradiated  sidebands  from  the  fundamental  laser  radia- 
tion. The  transmitted  intensity  of  the  sidebands  is  detected  with  high  sensitivity  in  a heterodyne 
receiver  in  which  a second  diplexer  combines  the  local  oscillator  and  sideband  beams  with  low 
loss  in  a corner-reflector-mounted  diode  mixer.  The  recovered  IF  is  amplified,  filtered,  and 
synchronously  detected.  The  new  sideband  system,  shown  in  Fig.  II-7 , is  significantly  more 
efficient  than  the  version  previously  discussed,  since  about  one-half  of  the  laser  radiation  is 
coupled  into  the  source  diode. 


Fig.  II— 7.  Tunable  sideband  submillimeter  spectrometer  witn  mgn-temperature  Mg  cell. 

Preliminary  measurements  with  the  system  have  included  the  observation  with  high  resolu- 
tion of  the  413  “*422  rotational  transition  in  D2O  at  692,245  ± 1 MHz  (692,240.5  ± 2.8  MHz  pre- 
dicted33 from  measurements  of  lower-frequency  transitions).  Also  observed  was  the  J = 5 — 6 

rotational  transition  in  CO  at  691,473  ± 1 MHz  (in  agreement  witn  the  measurements  of 

34  35 

Helminger,  De  Lucia,  and  Gordy  and  predictions  of  Kildal,  Eng,  and  Ross  based  upon  CO- 
laser  frequency  measurements  as  well  as  prior  data).  The  intensity  of  the  absorption  was  two- 
orders-of-magnitude  smaller  than  that  of  the  DzO  absorption.  The  local  oscillator  was  the 
432.6-ptm  line  of  HCOOH,  and  the  observations  were  made  with  a S/N  of  up  to  60  dB.  The  CO 
measurements  were  made  in  support  of  an  astronomical  program  to  detect  the  J = 5 — 6 transi- 
tion in  interstellar  dust  clouds. 
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The  high  S/N  attainable  makes  the  apparatus  attractive  for  the  study  of  still  weaker  transi- 
tion in  transient  species.  One  material  of  particular  interest  is  magnesium  in  its  metastable 
3P  state.  This  requires  more  sophisticated  sample  cells,  as  indicated  in  Fig.  II-7,  and  a spe- 
cialized high-temperature  cell  has  been  constructed  to  create  and  store  the  metastable  Mg.  Con- 
trolled small  quantities  of  Mg  effuse  from  an  externally  heated  source  chamber.  The  Mg  is  ex- 
3 

cited  to  the  P state  by  either  an  RF  or  DC  discharge.  The  quartz  windows  of  the  cell  are  cold 
and  protected  from  the  Mg  by  a buffer  gas.  Using  the  496-pm  line  of  CH,F,  we  have  generated 
a sideband  at  601.32  GHz  (corresponding  to  the  as  yet  unobserved  jPq  — JP^  fine-structure 
transition  in  Mg)  which  we  have  detected  with  a S/N  of  over  40  dB. 

We  are  increasing  the  sensitivity  of  the  spectrometer  by  amplitude  stabilization  of  the  sub- 
millimeter laser  and  by  using  Zeeman  modulation  of  the  Mg  cell  to  make  the  first  direct  obser- 
vation of  this  transition.  This  measurement  will  facilitate  the  locking  of  the  CH,F  laser  to  a 

36  * 

beam  of  metastable  Mg  for  use  as  a possible  submillimeter  frequency  standard. 

W.  A.  M.  Blumberg  D.  D.  Peck 
H.  R.  Fetterman  P.  F.  Goldsmith^ 
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III.  MATERIALS  RESEARCH 


A.  LASER  ANNEALING  OF  Se-IMPLANTED  GaAs 

Experiments  on  the  laser  heating  of  Si  and  GaAs  have  shown  that  amorphous  films  of  both 

materials  can  be  crystallized  by  this  technique  if  the  power  density  is  sufficient.1,2  In  addition, 

it  has  been  found  that  laser  heating  of  ion-implanted  Si  can  anneal  away  damage  and  produce 

good  electrical  activation  of  implanted  dopants,  in  some  cases  without  causing  appreciable 
3 

dopant  diffusion.  In  experiments  on  laser  annealing  of  ion-implanted  GaAs  layers,  backscat- 
tering  data  have  indicated  that  nearly  all  the  ion-implantation  damage  could  be  annealed  away 
and  that  the  dopant  species  were  in  substitutional  sites.4  However,  little  information  has  been 
published  on  the  electrical  characteristics  of  such  annealed  layers. 

We  have  now  carried  out  a study  of  the  effects  of  annealing  with  a CW  Nd:YAG  laser  on  the 
electrical  properties  of  Se-implanted  GaAs  layers.  Although  good  electrical  activation  has  been 
achieved,  we  have  observed  the  formation  of  (111)  defect  planes  which  degrade  the  electrical 
characteristics. 

The  samples  used  in  the  investigation  were  either  polished  slices  cut  from  semi-insulating 
Ci  -doped  crystals,  or  5-  to  10-pm-thick  epitaxial  layers  (S-doped  to  about  2 X 1013  cm-3)  on 
n+  substrates.  The  semi-insulating  samples  were  implanted  with  1 X 1014  cm-2  400-keV 
Se+  ions.  Control  samples  were  implanted  with  1 x 1014  cm-2  400-keV  Kr+  ions.  The  epi- 
taxial samples  were  in. planted  with  5 x 1012  cm-2  400-keV  Se+  ions.  Implant  temperatures 
were  25°,  250°,  and  300°C.  Some  of  the  samples  were  encapsulated  either  with  500  to  1000  A 
of  pyrolytic  Si^N4  deposited  prior  to  implantation,  or  with  SiOz  or  phosphosilicate  glass  (PSG) 
deposited  after  implantation. 


WHOLE  ASSEMBLY  SCANNED 
x-,y-  AND  l- DIRECTIONS 


Fig.III-1.  Schematic  representation  of  CW  Nd:YAG  laser  annealing  system. 

The  laser  annealing  system  is  shown  schematically  in  Fig.III-1.  The  multimode  output  of 
the  CW  Nd:YAG  laser  (operating  at  1.06  pm)  is  focused  by  the  cylindrical  lens  onto  the  sample. 
During  annealing,  the  sample  is  held  in  a forming  gas  (Ar/H2)  ambient  on  a platform  that  can 
be  resistively  heated  to  over  700°C.  The  platform  temperature  generally  used  was  580°C. 

This  temperature  is  high  enough  for  sufficient  bandgap  narrowing  to  allow  efficient  absorption 
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of  the  Nd:YAG  1.06-pm  radiation,  but  still  low  enough  so  that  no  thermal  dissociation  is  ob- 
served. The  platform  is  mounted  on  an  x,y,z  translational  stage  that  permits  the  sample  to  be 
scanned  under  the  laser  beam  in  the  focal  plane  and  normal  to  the  long  axis  of  the  slit  image. 

By  using  synchronous  motors,  scan  rates  from  0.05  to  19  mm/sec  have  been  obtained. 

Over  a hundred  implanted  GaAs  samples  have  been  laser-annealed  under  various  conditions. 
The  appearance  of  the  sample  surface  after  annealing  was  found  to  depend  on  the  laser  scan  rate 
and  power  density  P^.  For  each  scan  rate,  the  surface  appearance  changed  in  the  same  general 
manner  with  increasing  P^.  At  sufficiently  low  values  of  P^  the  entire  surface  remained  as 
smooth  as  before  annealing,  while  at  sufficiently  high  values  unencapsulated  surfaces  showed 
the  characteristic  features  produced  by  thermal  etching  due  to  the  loss  of  As  by  vaporization. 

The  threshold  P^  for  thermal  etching  increased  with  increasing  scan  rate,  from  less  than 
4 X 103  W/cm^  for  0.53  mm/sec  to  more  than  7 x 103  W/cm^  for  19  mm/sec.  For  each  scan 
rate,  over  an  intermediate  range  of  P^  values,  annealing  resulted  in  the  formation  of  character- 
istically oriented  straight  lines  on  the  otherwise  smooth  surface.  Such  features  have  not  been 
previously  reported  for  laser-annealed  samples  of  either  GaAs  or  Si.  The  values  of  PQ  at 
which  these  lines  were  first  observed  increased  with  increasing  scan  rate,  and  above  this  thresh- 
old the  density  of  lines  increased  with  increasing  PD.  The  highest  densities  of  lines  were  ob- 
served for  samples  annealed  at  the  highest  scan  rates.  Encapsulated  samples  generally  ex- 
hibited higher  densities  of  lines  than  unencapsulated  samples  annealed  under  the  same  conditions, 
especially  at  high  scan  rates. 


Fig.III-2.  Micrographs  of  top  surface  [2°  off  (100)  plane  in  (110)  direction] 
and  a cleaved  (011)  cross  section  of  a laser-annealed  GaAs  sample  showing 
lines  due  to  defects  produced  by  annealing.  Schematic  diagram  locates  lines 
three -dimensionally  on  sample.  Dark  arrow  in  diagram  indicates  laser  scan 
direction.  Defects  have  tentatively  been  identified  as  (111)  slip  planes. 

The  characteristics  of  the  surface  lines  are  illustrated  by  Fig.  III-2,  which  shows  two  optical 

micrographs  of  an  epitaxial  sample  annealed  at  the  highest  scan  rate  (19  mm/sec)  using  Pn  = 

3 2 ^ 

7.3  X 10  W/cm  . The  upper  micrograph  shows  the  top  surface  of  the  sample,  which  was  ori- 
ented 2*  off  the  (100)  plane,  while  the  lower  one  shows  an  (011)  plane  cleaved  perpendiculu'-  to 
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the  top  surface  and  stained  with  a potassium  ferricyanide  solution.  The  two  micrographs  have 
been  positioned  to  align  the  corresponding  features  on  the  top  surface  and  cleaved  cross  section. 
Figure  III-2  also  includes  a schematic  diagram  locating  these  features  three -dimensionally  with 
respect  to  the  sample,  which  was  formed  by  the  epitaxial  growth  of  an  n+  buffer  layer  and  an 
upper  n-layer  on  an  n+  GaAs  substrate. 

The  upper  micrograph  of  Fig.  III-2  shows  two  pairs  of  parallel  straight  lines  ("tracks")  on 
the  top  surface  of  the  sample  - one  pair  at  an  angle  of  1 • to  the  cleaved  edge  at  the  intersection 
between  the  (100)  and  (Oil)  planes,  the  other  pair  also  at  a 1*  angle  to  this  edge  but  in  the  op- 
posite sense.  The  angle  of  2°  between  the  two  pairs  of  tracks  results  from  the  2°  deviation  of 
the  top  surface  from  (100)  toward  (110).  As  can  be  seen  in  the  micrograph  of  the  (Oil)  cleaved 
face,  the  tracks  propagate  from  the  surface  through  both  epitaxial  layers  into  the  substrate. 

On  the  basis  of  their  orientation,  we  have  tentatively  identified  the  surface  tracks  as  intersec- 
tions between  the  (100)  surface  and  two  sets  of  (111)  slip  planes,  which  intersect  with  each  other 
at  70°  in  the  plane  of  the  cleaved  (Oil)  face.  Since  the  slip  planes  propagate  from  the  surface 
at  an  angle  of  55°,  their  density  is  always  highest  at  the  surface.  Except  for  the  tracks,  the 
top  surface  is  as  smooth  as  it  was  before  annealing,  and  the  interface  between  the  upper  epi- 
taxial layer  and  the  buffer  layer  is  also  unchanged,  indicating  that  no  melting  occurred. 

Formation  of  the  slip  planes  is  presumably  caused  by  thermal  stress  that  occurs  during 
laser  annealing.  This  mechanism  is  consistent  with  the  observation  that  the  threshold  PD  for 
producing  the  surface  lines  increased  with  increasing  laser  scan  rate,  since  the  Pp  required 
to  achieve  a given  surface  temperature  also  increases  with  scan  rate. 

The  presence  of  slip  planes  has  a strong  adverse  effect  on  the  electrical  properties  of  laser- 
annealed  GaAs.  Tli is  effect  can  be  seen  from  the  data  listed  in  Table  IH-1,  which  gives  the  ef- 
fective sheet  carrier  concentration  Ns  and  effective  mobility  p for  representative  samples  pre- 
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pared  by  implantation  of  1 x 10  cm  Se  ions  directly  into  semi-insulating  substrates  and 

then  annealed  under  various  conditions.  Significant  trends  are  illustrated  by  Fig.  III-3,  in  which 

Ng  is  plotted  against  PD  for  a number  of  unencapsulated  samples  that  were  Se-implanted  at 

300°C,  then  subjected  to  single  laser  scans  at  four  different  scan  rates.  The  values  of  N and 

s 

u for  control  samples  that  were  Se-implanted  at  300 °C  and  then  heated  to  580 °C  but  not  laser 
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annealed  were  9 x 10  cm  and  980  cm  /V-sec,  respectively.  For  such  control  samples 
implanted  at  room  temperature,  the  resistivity  was  too  high  to  permit  Hall  measurements.  All 
the  Kr  -implanted  control  samples  remained  semi-insulating  following  laser  annealing. 

The  data  of  Fig.  HI— 3,  which  are  restricted  to  samples  that  did  not  have  significant  thermal 
etching,  show  an  overall  decrease  in  Ng  with  increasing  scan  rate.  At  the  highest  scan  rate, 
over  the  range  of  PD  values  we  examined,  Ng  actually  decreases  with  increasing  PD.  On  the 
other  hand,  at  the  lower  scan  rates  Ng  increases  with  increasing  P^.  These  results  may  be 
explained  by  assuming  that  laser  annealing  can  have  two  opposing  effects:  removing  ion- 
implantation  damage  and  producing  slip  planes.  Therefore,  it  can  increase  Ng  by  removing 
the  compensating  defects  due  to  implantation  damage,  but  it  can  also  decrease  Ng  by  introduc- 
ing compensating  defects  due  to  slip  planes.  For  the  two  lower  scan  rates  shown  in  Fig.  III-3, 

Ns  increases  with  increasing  P^  because  the  increased  removal  of  implantation  damage  out- 
weighs the  rather  limited  formation  of  slip  planes  in  these  samples.  The  sample  annealed  at 
0.53  mm/sec  with  P^  = 3.0  x 10^  W/cm^  has  the  highest  Ng  that  we  have  obtained  for  a singly 
scanned  sample,  Ng  = 2.7  X lO*'*  cm  ^ with  pg  = 2100  cm^/V-sec.  By  comparison,  similarly 
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Fig.  III-3.  Effective  sheet  carrier 
concentration  vs  laser  power  density 
at  scan  rates  from  0.53  to  19  mm/sec 
for  GaAs  samples  implanted  with  1 X 

1014  cm  2 400-keV  Se+  ions  (implant 
temperature  300 °C). 


implanted  samples  that  are  thermally  annealed  at  900°C  for  15  min.  with  a pyrolytic  Si,N.  en- 
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capsulant  typically  have  Ng  = 5.0  x 10  cm  and  = 2300  cm  /V-sec.  Further  improve- 
ment can  be  expected  at  0.53  mm/sec  by  increasing  P^,  provided  that  the  threshold  for  thermal 
etching  is  not  significantly  exceeded. 

For  the  samples  of  Fig.  III-3  that  were  annealed  at  a scan  rate  of  9.5  mm/sec,  all  of  which 
exhibited  significant  densities  of  surface  lines,  Ng  was  essentially  independent  of  PD  over  the 
range  used,  indicating  that  the  two  opposing  effects  of  laser  heating  canceled  each  other.  Of 
the  samples  annealed  at  19  mm/sec,  the  two  annealed  at  the  lower  values  of  P^  did  not  exhibit 
surface  lines  but  both  had  the  same  low  value  of  Ng.  Probably  the  surface  temperatures 
achieved  at  those  power  levels  were  not  sufficient  to  remove  very  much  of  the  implantation 
damage.  For  the  third  sample,  which  was  annealed  at  6 x 10  W/cm  , not  many  surface  lines 
were  observed,  but  these  apparently  outweighed  any  increased  removal  of  implantation  damage. 
The  use  of  still  higher  power  levels  resulted  in  a rapid  increase  in  the  density  of  the  surface 
lines,  and  the  results  of  van  der  Pauw  measurements  became  too  anisotropic  to  be  meaningful. 

The  anisotropy  in  the  van  der  Pauw  results  is  another  consequence  of  the  formation  of  slip 
planes.  This  is  demonstrated  by  the  results  shown  in  Fig.  III-4,  which  includes  two  I-V  curves 
obtained  for  different  electrode  pairs  in  measurements  on  a cloverleaf  mesa  etched  on  the  sur- 
face of  a sample  that  had  been  laser-annealed  under  conditions  that  produced  a high  density  of 
slip  planes.  The  schematic  diagram  of  the  cloverleaf  in  Fig.  Ill -4  shows  the  configuration  of 
the  surface  lines  in  relation  to  the  electrodes  (only  three  of  the  lines  are  indicated).  The  I-V 
curve  designated  as  (a),  which  was  obtained  for  electrodes  i and  2 or  3 and  4,  is  linear  and 
indicates  low  resistivity.  Curve  (b),  for  electrodes  1 and  4 or  2 and  3,  is  characteristic  of 
back-to-back  diodes.  In  the  first  case,  the  current  paths  between  the  electrodes  do  not  inter- 
sect the  slip  planes,  while  in  the  second  case  they  do.  The  characteristics  of  curve  (a)  thus 
show  that  the  laser  annealing  removed  the  ion  implantation  damage  in  the  regions  between  the 
planes,  while  those  of  curve  (b)  show  that  the  planes  act  as  barriers  in  n-type  GaAs. 

The  effects  of  the  slip  planes  on  the  electrical  properties  of  laser-annealed  samples  are 
further  shown  by  C-V  carrier-concentration  profiles  on  several  ion-implanted  GaAs  epitaxial 
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Fig.  III-4.  Current-voltage  characteristics  between  different  contacts  to  an  etched 
mesa  cloverleaf  on  laser-annealed  GaAs  with  a high  density  of  defect  lines.  Insert 
schematically  shows  etched  cloverleaf  with  three  such  lines,  (a)  I-V  character- 
istics parallel  to  these  lines,  i.e.,  between  contacts  1 and  2 or  3 and  4;  (b)  across 
the  lines,  i.e.,  between  contacts  1 and  4 or  2 and  3. 


Fig.III-5.  Effective  carrier  concentration 
vs  depth  obtained  from  C-V  measurements 

on  GaAs  samples  implanted  with  5 x 1012 

cm  2 400-keV  Se+  ions  (implant  temper- 
ature 300  °C).  Solid  line  represents  data 
for  an  encapsulated  sample  that  was  ther- 
mally annealed  at900°C  for  15  min.  Other 
three  samples  were  mounted  on  substrate 
platform  which  was  heated  to  580°C. 


layers.  These  layers  had  an  initial  carrier  concentration  of  2 X lO1^  cm  3 and  were  implanted 
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with  5 x 10  cm  Se  ions  at  300  °C.  Figure  III-5  shows  the  effective  carrier -concent ration 
(nefj)  profiles  for  four  such  samples.  The  solid  line  represents  data  obtained  for  a sample  that 
was  encapsulated  with  SijN4  and  thermally  annealed  at  900°C  for  15  min.  The  other  three  sam- 
ples were  not  encapsulated.  In  separate  experiments  using  the  laser-annealing  system,  each 
was  mounted  on  the  sample  platform  which  was  heated  to  580 °C.  Hie  sample  not  subjected  to 
laser  heating  exhibited  no  activation.  For  the  sample  subjected  to  a single  laser  scan  at 
9.5  mm/sec  with  = 4.9  X 10  W/cm  , slip  planes  were  observed  but  activation  occurred. 
However,  the  peak  ne^  was  much  lower  than  that  obtained  for  the  thermally  annealed  sample. 

In  addition,  n .»  in  the  unimplanted  region  of  the  epitaxial  layer  decreased  from  its  original 
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level  of  about  2 x 10  cm”  to  about  10  3 cm”  , presumably  because  the  slip  planes  prop- 
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agated  into  this  region.  For  the  sample  annealed  at  0.05  mm/sec  with  Pp  = 3 X 10  W/cm  , 
the  peak  neff  was  higher  than  in  the  other  laser-annealed  sample,  but  still  lower  than  the  value 
obtained  by  thermal  annealing.  In  this  sample,  the  slip  planes  probably  account  for  the  deep 
minimum  in  n^.  but  they  apparently  propagated  to  a depth  of  only  about  4 pm,  since  ne^  re- 
turned to  its  initial  value  at  about  this  depth. 

Reflection  high-energy-electron  diffraction  (RHEED)  indicates  that  the  surfaces  of  samples 
implanted  at  250°  or  300°C  are  single  crystal  both  before  and  after  laser  annealing.  Samples 
implanted  at  room  temperature  were  found  to  be  amorphous  prior  to  annealing,  and  single  crys- 
tal afterwards.  It  was  further  noted  that  these  RHEED  observations  show  good  single  crystal- 
linity even  when  a large  number  of  slip  planes  are  present.  This  result  indicates  that  in  order 
to  characterize  laser-annealed  samples,  crystallographic -type  observations  should  be  supple- 
mented by  electrical  measurements. 

In  summary,  one  laser  scan  is  sufficient  to  obtain  fairly  good  activation  of  GaAs  samples 

implanted  at  250°  or  300°C.  The  best  results  obtained  were  for  low  scan  rates.  At  0.53  mm/ 
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sec,  the  best  results  for  a single  scan  were  N = 2.7  x 10  cm  and  p = 2100  cm  /V-sec. 
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On  samples  scanned  ten  times  at  0.53  mm/sec,  the  best  measured  N was  5.9  X 10  cm  . 

2 ® 

The  corresponding  sheet  mobility  was  only  1300  cm  /V-sec,  however,  and  is  most  likely  lim- 
ited by  slip-plane  formation.  For  fast  scan  rates,  multiple  scans  usually  caused  more  defect 
planes,  and  Hall  measurements  became  quite  anisotropic. 

On  samples  implanted  at  room  temperature,  one  scan  did  not  generally  result  in  appreciable 

activation,  but  ten  scans  did.  The  best  results  obtained  for  a sample  implanted  at  room  tempera- 
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ture  and  scanned  ten  times  at  0.53  mm/sec  were  Ng  = 2.2  x 10  cm  and  pg  = 2200  cm  /V-sec. 

Overall,  the  results  suggest  that  very  good  electrical  characteristics  may  be  possible  if  the 
formation  of  slip  planes  can  be  prevented  or  minimized.  Since  these  planes  are  probably  caused 
by  thermal  stress  encountered  in  laser  heating,  annealing  procedures  such  as  lower  scan  rates 
and/or  higher  substrate  platform  temperatures,  which  reduce  thermal  stress,  should  lead  to 
some  improvement.  It  may  also  be  possible  that  slip-plane  formation  can  be  reduced  if  surface 
orientations  other  than  (100)  are  used.  There  is  probably  no  obvious  advantage  to  using  pulsed- 
laser  heating,  since  large  thermal  stresses  are  incurred  in  this  mode  of  operation. 

J.  C.C.  Fan  C.O.  Bozler 

J.  P.  Donnelly  R.  L.  Chapman 

B.  ACTIVATION  ENERGY  OF  Cr  ACCEPTORS  IN  InP 

Infrared  diode  lasers  and  avalanche  photodiodes  utilizing  GalnAsP  epitaxial  layers  grown 
on  InP  substrate  are  promising  sources  and  detectors  for  fiber  optics  communications  at  1.1  to 


Fig.  III-6.  Dependence  of  (left 

ordinate)  and  resistivity  p (right  ordinate) 
on  reciprocal  absolute  temperature  l/T  for 
an  InP sample  from  a Cr-dopedLEC  boule. 


Fig.  III-7.  Calculated  dependence  of  carrier 
concentration  in  InP  at  295  K,  n,Q,,  on  N.  / 


the  ratio  of  deep  acceptor  (Fe  or  Cr) 

concentration  to  difference  between  total  do- 
nor concentration  and  concentration  of  ac- 
ceptors other  than  Fe  or  Cr. 
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1.3  pm.  As  previously  reported,5  we  have  been  growing  InP  single  crystals  by  the  liquid- 
encapsulated  Czochralski  (LEC)  method  to  provide  substrate  material  for  use  in  the  develop- 
ment of  these  devices.  Semi-insulating  crystals  with  resistivities  of  ~108  n-cm  have  been  ob- 
tained by  doping  with  Fe,  a deep  acceptor.  By  measuring  the  Hall  coefficient  (R^)  as  a function 
of  temperature  for  a sample  from  one  of  these  crystals,  the  activation  energy  for  thermally 
exciting  electrons  from  the  Fe  acceptor  level  to  the  conduction  band  was  found  to  be  0.65  eV 
(see  Ref.  5). 

Like  Fe,  Cr  is  known  to  be  a deep  acceptor  in  InP,  since  resistivities  of  103  to  105  Si-cm 
have  been  obtained  by  Cr  doping  of  material  that  would  otherwise  be  n-type  with  low  resistivity.6 
However,  these  values  are  much  lower  than  those  achieved  by  Fe  doping.  In  order  to  establish 
the  reason  for  this  difference,  we  have  determined  the  activation  energy  for  electron  excitation 
from  the  Cr  acceptor  level.  As  in  the  case  of  Fe,  this  determination  was  made  by  measuring 
Rjj  as  a function  of  temperature.  The  measurements  were  made  on  a Cr-doped  sample,  cut 
from  a crystal  grown  by  the  LEC  method,  with  room-temperature  resistivity  (p)  of  5.3  X 
10  fi-cm,  carrier  concentration  (n  = l/RHe)  of  2.6  X 1010  cm-3,  and  Hall  mobility  (p  = RH/p) 
of  4600  cm2/V-sec.  The  results  are  shown  in  Fig.  III-6,  where  RHT3/2  and  p are  plotted  on 
logarithmic  scales  vs  l/T.  The  activation  energy  obtained  from  the  R^T3^2  plot  is  0.39  ± 

0.01  eV,  in  good  agreement  with  the  value  of  0.37  eV  obtained  by  Pande  and  Roberts7  as  the 
thermal  activation  energy  for  the  flow  of  current  through  a Cr-doped  sample  in  the  ohmic  cur- 
rent range.  For  the  space-charge-limited  range,  these  authors  found  an  activation  energy  of 

0.20  eV.  They  explained  their  observations  in  terms  of  a deep  donor,  deep  acceptor  model  with 
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a total  deep  acceptor  concentration  of  <10  cm  . This  model  is  inconsistent  with  our  results, 
since  it  implies  a residual  net  donor  concentration  that  is  also  <1012  cm’3,  while  we  have  never 
observed  values  of  n30Q  below  about  5 x 1014  cm’3  for  nominally  undoped  samples  of  InP. 

The  molar  Cr  concentration  in  the  charge  used  for  our  Cr-doped  boule  was  4500  ppm.  As- 
suming that  Cr  is  a simple  deep  acceptor  with  an  activation  energy  of  0.39  eV  and  that  the  re- 
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sidual  donor  concentration  is  3 X 10  cm  , from  the  measured  value  of  n30Q  for  a first -to- 

freeze  sample  we  have  calculated  the  concentration  of  Cr  incorporated  in  the  sample  to  be 
2.7  ppm,  corresponding  to  a distribution  coefficient  of  6 X 10"4.  This  value  falls  within  the 
range  of  1 to  6 X 10  found  by  means  of  mass  spectrographic  analysis. 

For  wide-bandgap  semiconductors  containing  excess  residual  donors  (i.e.,  > 0,  where 

Np  is  the  difference  between  the  total  donor  concentration  and  the  concentration  of  shallow  ac- 
ceptors), very  high  resistivities  can  be  obtained  by  incorporating  acceptor  impurities  with  levels 
located  near  the  center  of  the  bandgap  because  the  Fermi  level  is  shifted  from  the  vicinity  of  the 
conduction  band  to  the  vicinity  of  the  deep  levels,  greatly  reducing  the  carrier  concentration. 

The  shift  of  the  Fermi  level  increases  with  increasing  separation  between  the  deep  levels  and 
the  conduction  band  Ec>  and  also  as  the  ratio  of  the  deep  acceptor  concentration  N . to  N*  in- 

Ad  D 

creases.  For  InP  doped  with  Fe  or  Cr,  these  trends  are  shown  in  Fig.  III-7  where  calculated 

values  of  n at  295  K are  plotted  against  /N|j.  For  doping  with  Fe  (acceptor  level  0.65  eV 

d 

below  Ec),  n decreases  precipitously  as  soon  as  the  Fe  concentration  exceeds  N*j,  reaching 

values  of  10  and  10  cm  for  /N^  ratios  of  1.2  and  2.6,  respectively.  In  contrast,  for  dop- 

d 

ing  with  Cr  (acceptor  level  0.39  eV  below  E ),  n does  not  reach  10 10  cm'3  until  N /n*  = 50 

8 3 Ad  D ’ 

and  is  still  9 x 10  cm  even  for  /N*j  = 500.  Thus,  the  highest  resistivities  obtained  for 

d 

Cr-doped  InP  are  much  lower  than  those  achieved  by  Fe  doping  because  the  Cr  acceptor  level 


is  considerably  closer  to  the  conduction  band  than  is  the  Fe  level,  not  because  Cr  is  less  soluble 
in  InP  than  Fe.  Even  if  the  solubility  of  Cr  were  the  same  as  that  of  Fe,  which  is  reported  to 
be  2.5  x 10  cm  3 (Ref.  9),  and  N|j  were  only  1 x 1015  cm'3,  n would  still  be  ~2  x 109  cm'3 
at  the  Cr  solubility  limit,  compared  with  the  values  of  ~107  cm'3  we  have  obtained  for  Fe-doped 
samples. 

G.  W.  Iseler 
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IV.  MICROELECTRONICS 


A.  INTEGRATED  MIXER  MODULE  FOR  SUBMILLIMETER-WAVELENGTH  OPERATION 

A GaAs  monolithic  integrated  circuit  capable  of  receiving  radiatively  coupled  submillimeter  - 
wavelength  radiation  has  been  fabricated.  Figure  rv-i  is  a schematic  of  the  circuit  elements 
which  have  been  integrated  on  this  module.  The  radiation  is  received  by  a slot  antenna  and 
coupled  to  a surface -oriented  mixer  diode ^ by  an  appropriate  section  of  coplanar  transmission 
line.  An  integrated  bypass  capacitor  completes  the  mixer  circuit  by  providing  a low  impedance 
to  radio  frequencies  and  a high  impedance  to  intermediate  frequencies.  The  realization  of  the 
circuit  of  Fig.IV-1  is  shown  in  Fig.IV-2(a).  The  design  of  the  planar  slot  antenna  is  complicated 
by  the  presence  of  the  high-resistivity  GaAs.  Free-space  radiation  patterns,  even  from  a simple 
metal  antenna  fabricated  on  the  high-resistivity  GaAs,  show  complex  lobe  structure  unacceptable 
for  most  applications.  Therefore,  the  module  of  Fig.rv-2(a)  has  been  designed  so  that  the  radi- 
ation is  coupled  through  the  high-resistivity  GaAs  substrate.  The  radiation  pattern  obtained  at 
350  GHz  is  shown  in  Fig.IV-2(b).  It  is  nearly  identical  to  that  obtained  from  a frequency -scaled 
model  at  5 GHz . The  pattern  was  obtained  in  both  cases  by  measuring  the  video  response  of  the 
surface -oriented  diode  as  the  direction  of  the  incident  radiation  was  changed.  The  theoretically 
expected  radiation  pattern  from  a half -wavelength  dipole  is  also  shown  in  Fig.  IV -2(b)  for 
comparison . 
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LINE  DIODE  CAPACITOR 

Fig.IV-1.  Circuit  diagram  of  integrated 
mixer  module  designed  for  quasi-optical 
operation  at  350  GHz. 


Fig.IV-2.  Millimeter -wave  mixer  mod- 
ule. (a)  Configuration  of  device  and  cir- 
cuit elements  on  module.  Radiation  is 
coupled  through  GaAs  as  shown,  (b)  Ra- 
diation patterns  at  350  GHz. 
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Though  the  measured  radiation  pattern  does  contain  some  interference  from  the  coplanar 
transmission-line  feed,  it  allows  the  near-millimeter-wavelength  radiation  to  be  radiatively 
coupled  to  the  device.  The  conversion  loss  measured  at  350  GHz  with  this  mixer  was  approxi- 
mately 20  dB,  as  compared  with  a 7 .5-dB  loss  obtained  with  the  frequency-scaled  model. 
Network -analyzer  measurements  have  indicated  that  the  comparatively  poor  performance  of  the 
mixer  module  at  350  GHz  is  due  to  loss  in  the  RF  bypass  capacitors.  A mixer  module  which 
incorporates  improved  bypass  capacitors  is  currently  being  fabricated. 

R.A.  Murphy 
G.D.  Alley 
B.  J.  Clifton 

B.  CHARGE-COUPLED  DEVICES:  IMAGERS 

The  100-  x 400 -element  CCD  imagers  fabricated  for  the  GEODSS  (Ground  Electro-Optical 
Deep  Space  Surveillance)  program  are  required  to  have  as  high  a sensitivity  as  possible.  This 
implies  both  high  quantum  efficiency  and  low-noise  charge  transfer  and  charge  detection.  In 
this  report,  we  describe  some  experimental  results  on  noise  measurements  related  to  charge 
detection  and  signal  processing.  Using  correlated  double  sampling  (CDS),  we  have  achieved 
noise-equivalent  signals  from  the  detection  circuitry  of  10  electrons  at  a 400 -kHz  data  rate. 


Fig.IV-3.  Diagram  of  system  used  in  noise  measurements. 


A diagram  of  the  system  used  to  evaluate  CCD  noise  is  shown  in  Fig.IV-3.  Charge  detec- 
tion on  the  chip  is  done  by  a conventional  floating-diffusion  output  circuit  consisting  of  a reset 
transistor  T4  and  output  transistor  T2-  The  signal  for  T2  is  buffered  by  the  emitter  follower  T3 
which  is  a discrete  2N930  transistor  located  in  the  device  package  next  to  the  CCD  chip2  along 
with  the  load  resistors  for  T2  and  Tj.  The  signal  is  then  amplified  by  a low-noise  amplifier 
having  an  input  noise  voltage  of  less  than  2 nV  Hz-1/2.  The  gain  A is  adjusted  to  a level, 
typically  50,  where  noise  contributions  from  the  remaining  circuits  are  insignificant.  Spectral 
measurements  of  the  noise  from  T2  and  T3  are  made  at  this  point.  The  signal -processing  chain 
continues  with  a low-pass  filter  formed  by  Rp  and  Cp  to  band-limit  the  noise,  followed  by  a 


circuit  which  performs  correlated  double  sampling.  “ The  mean  and  variance  of  the  sampled 
voltage  signal  are  measured  by  first  digitizing  a sequence  of  signal  samples  (typically  1024)  to 
10-bit  resolution  in  a waveform  digitizer^  and  transmitting  the  data  to  a programmable  calcu- 
latort  for  computation.  The  calculator  can  plot  on  a display  scope  the  distribution  function  of 
the  voltage  samples  together  with  a comparison  plot  of  a Gaussian  curve  of  the  same  mean  and 
standard  deviation  to  verify  the  Gaussian  nature  of  the  noise. 

The  noise -equivalent  signal  in  electrons  is  then  determined  by  the  rms  noise  voltage  after 
signal  processing  divided  by  the  voltage/electron  conversion  ratio.  The  term  usually  applied  to 
this  ratio  is  responsivity  R expressed  in  volts/electron,  and  is  given  by 


where  G is  the  voltage  gain  between  the  gate  of  T2  and  the  point  in  the  circuit  where  the  signal 
is  measured,  and  CQ  is  the  capacitance  at  the  output-diode  node.  Low-noise  operation  then 
dictates  low  CQ  and  low  noise  from  T2 . By  careful  design  of  the  output  circuit,  the  capacitance 
CQ  has  averaged  0.050  pF  or  a factor-of-4  lower  than  values  typical  of  CCDs  described  in  the 
literature.  Corresponding  responsivities  measured  at  the  emitter  of  Tj,  where  G = 0.75,  have 
been  in  the  range  of  2.3  to  2.5  pV/electron. 


DEVICE  NO.  118-14 


Fig. IV -4.  Measurement  of  noise  spectral 
voltage  of  output  MOSFET  and  emitter  fol- 
lower. MOSFET  is  operating  in  buried - 
channel  mode. 
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Under  ideal  circumstances,  the  off-chip  signal  processing  should  result  in  noise  levels 

4 5 

determined  solely  by  the  thermal  noise  of  T2-  It  has  been  shown  by  Brodersen  and  Emmons  ' 
that  the  noise  of  T2  can  be  significantly  reduced  if  this  transistor  is  subjected  to  the  same  buried - 
channel  implant  as  the  remainder  of  the  device,  and  if  it  is  operated  in  a buried -channel  mode 
(i.e.,  the  channel  current  flow  is  kept  away  from  the  oxide-silicon  interface).  The  normally  high 
l/f  noise  is  significantly  reduced  in  this  case.  In  our  process,  T2  receives  the  buried-channel 
implant,  and  we  have  observed  greatly  reduced  noise  output  when  operating  T2  buried  channel. 

A typical  measurement  of  noise  spectral  voltage  for  the  T2-T,  combination  is  shown  in  Fig. IV -4, 


t Biomation  Model  1010. 
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and  it  is  seen  that  significant  l/f  noise  still  occurs  at  frequencies  below  50  kHz.  The  measured 
slope  of  the  low-frequency  noise  voltage  exhibits  a l/f^'3^  dependence  for  this  device,  in  con- 
trast to  the  usual  l/f1/2  relationship.  The  high-frequency  noise  level  of  12  n\/*J~Hz  is  somewhat 
higher  than  the  theoretical  value  of  8.6  nV/V  Hz. 

The  signal  processing  performed  by  the  low-pass  filter  and  correlated  double  sampling  is 
needed  to  fully  realize  the  low-noise  capability  of  this  system.  The  low-pass  filter  removes  high- 
frequency  noise  which  would  be  aliased  into  the  sampled  data  passband  below  fc/2,  where  fc  is 
the  clock  frequency.  The  optimum  corner  frequency  f.  of  this  filter  has  been  shown  to  be4<5 

Li 

fL  = 4V3  ' (IV-2) 

The  video  clamping  performed  by  coupling  capacitor  C and  switch  T4  in  Fig.IV-3  removes  reset 
noise  from  the  signal  when  the  reset  gate  of  the  CCD  is  being  clocked,  and,  in  addition,  helps 
suppress  low-frequency  noise  such  as  l/f  noise.  Reset  noise  is  given  by  (kTC  J1/2 /q,  and  would 

2 o 

be  90  electrons  in  our  case.  The  mean-square  output -noise  voltage  VQ  from  the  low-pass  filter/ 
CDS  processor  in  terms  of  the  input-noise  spectral  voltage  V?  is5 

V2  = 2 [R(t>  - R(o)] 

where 

v? 

R(t)  = \ 5^-5- (cos  2irff)  df  (IV-3) 

Jo  1 + f2/f 2 

is  the  autocorrelation  function  of  the  noise,  and  r is  the  time  between  clamp  and  sample  pulses 
and  is  typically  1 / (2f c> . 

The  circuitry  described  above  permitted  the  observation  of  very  low  noise  levels.  As  a 
check  on  the  accuracy  of  the  noise  measurements,  the  CCD  was  illuminated  with  an  incandescent 
bulb  and  the  resulting  signal  and  noise  levels  were  measured  at  the  output  of  the  CDS  circuit. 

For  this  measurement,  only  the  405-stage  output  register  of  the  imager  was  clocked  to  minimize 
dark  current.  In  theory,  optically  generated  charge  packets  having  an  average  of  n electrons/ 
packet  would  have  a photon  shot  noise  level  of  (n)1/2  electrons.  The  experimental  results  in 
Fig.IV-5  show  agreement  between  theory  and  experiment  to  within  5 percent. 


Fig. IV -5.  Measured  noise  level  on  opti- 
cally generated  signal.  Theoretical  pho- 
ton shot  noise  level  in  electrons  is  equal 
to  square  root  of  average  signal  level  and 
is  shown  as  solid  line. 
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Fig. IV -6.  Measured  amplifier  noise  from  CCD  with  no  sig- 
nal charge.  Theoretical  curve  is  based  on  measured  spectral 
noise  in  Fig.IV-4  and  includes  a calculation  assuming  the  data 
in  Fig.IV-4  had  no  i/f  component. 

Measurements  were  made  of  the  system  noise  with  the  CCD  clocked  in  the  normal  fashion, 
but  with  the  output  gate  biased  off.  In  this  case,  dark-current  charge  is  prevented  from  being 
delivered  to  the  CCD  output  circuit,  and  the  system  noise  should  be  only  that  of  the  output  cir- 
cuit. Measurements  made  at  clock  frequencies  from  10  kHz  to  1 MHz  are  shown  in  Fig.IV-6. 

At  each  frequency,  the  low-pass  filter  was  adjusted  according  to  Eq.(IV-2).  The  values  range 
from  about  7 electrons  at  10  kHz  to  20  electrons  at  1 MHz,  and  represent  the  lowest  reported 
noise -equivalent  signals  from  a floating  diffusion  output  structure.  Plotted  in  the  same  figure 
is  a theoretical  curve  using  Eq.(IV-3)  and  the  measured  noise  data  from  Fig.IV-4.  In  addition, 
the  calculation  was  also  performed  (dashed  curve)  assuming  the  measured  noise  data  had  only 
the  high-frequency  white  noise  and  no  l/f  component.  This  curve  suggests  that  noise  levels 
below  one  electron  may  be  possible  at  low  clock  frequencies  with  carefully  made  low-noise  out- 
put MOSFETs.  The  measured  data  are  about  a factor-of-2  higher  than  theory.  The  major  dif- 
ference in  the  measurement  conditions  of  Figs.IV-4  and  IV -6  is  that  in  the  former  the  CCD  out- 
put is  DC -biased  to  its  operating  point  with  all  clocks  off,  while  in  the  latter  the  CCD  is  clocked. 
Earlier  measurements  of  noise  were  much  higher,  and  the  problem  was  traced  to  a commercial 
pulse  generator  used  for  the  reset  gate  which  had  fairly  high  levels  of  noise  (several  thousand 
nV/VHz).  Generators  having  much  lower  noise  levels  were  then  used,  and  much  lower  CCD 
noise  levels  were  achieved.  The  remaining  discrepancy  between  theory  and  experiment  may 
still  be  due  to  pulse-generator  noise. 

B.E.  Burke 
W.H.  McGonagle 

C.  CHARGE-COUPLED  DEVICES:  SAW/CCD  BUFFER  MEMORY 

The  SAW/CCD  buffer  memory  is  described  in  detail  in  earlier  reports. ’7  The  device  is 
comprised  of  a set  of  sampling  fingers  on  a CCD  die  which  is  in  close  proximity  (typically  300  nm) 
to  a SAW  delay  line.  The  piezoelectric  RF  fields  associated  with  a traveling  surface  wave  are 
sampled  and  held  by  the  fingers,  and  the  resulting  charge  pattern  is  parallel -transferred  to  a 
CCD  shift  register.  The  overall  frequency  response  of  the  SAW/CCD  buffer  memory  is  affected 
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by  the  frequency  response  of  the  input  transducer  and  the  fact  that  the  sampling  fingers  have 
finite  width.  The  frequency  response  of  the  input  transducer  is  well  understood8  and  will  not  be 
discussed  here.  In  this  report  we  present  a theoretical  analysis  of  the  effect  of  finger  width  on 
frequency  response  and  compare  the  results  of  this  analysis  with  the  measured  response  of  the 
prototype  device. 

The  finite  width  of  the  sampling  fingers  causes  the  frequency  response  to  roll  off  gradually 
at  higher  frequencies  within  the  input  frequency  band.  If  we  assume  that  each  finger  samples 
the  average  of  the  sinusoidally  varying  piezoelectric  field  across  its  width,  then  this  component 
of  the  frequency  response  is  given  by  the  integral  of  the  input  waveform  over  the  width  of  the 


sampling  finger  divided  by  the  finger  width.  For  a sinusoidal  input,  this  integral  is  identical  to 
the  Fourier  transform  of  a rectangular  pulse?  For  a finger  of  width  W,  the  transfer  function 


is  given  by 


™ „ sin  (irW/X) 
T-F-  - irW/x 


(IV -4) 


Here,  X is  the  wavelength  (or  the  inverse  of  the  spatial  frequency)  of  the  piezoelectric  field  and 
is  equal  to  V /f,  where  V is  the  velocity  of  the  surface  acoustic  wave  and  f is  the  input  radio 


frequency.  This  transfer  function  is  more  useful  if  it  is  expressed  in  terms  of  f,  fg,  and  SQ: 


T.F.  = 


sin  ir(W/SJ  (f/f  ) 


*(W/Sc)  (f/fg) 


(IV— 5 ) 


where  Sc  is  the  spatial  period  of  the  sampling  fingers,  and  fg  = Va/Sc  is  an  effective  sampling 
frequency  discussed  in  an  earlier  report.6  A computer  program  has  been  written  which  calcu- 
lates the  theoretical  frequency  response  of  the  SAW/CCD  buffer  memory  by  combining  the  above 
equation  with  the  theoretical  response  of  the  input  transducer  and  its  matching  network. 


Fig.IV-7.  Frequency  response  of 
buffer  memory. 


The  frequency  response  of  the  first  prototype  device  is  shown  in  Fig.IV-7.  The  measured 
response  falls  off  inside  the  pass  band  of  the  input  transducer  as  predicted  by  theory.  (The 
artifact  in  the  experimental  data  around  98  MHz  is  the  result  of  measurement  errors,  which  can 
be  as  high  as  1 dB.)  The  theoretical  curve  in  the  figure  is  calculated  for  W = 18  pm,  which  is 
equal  to  the  measureo  finger  width  of  16  pm  plus  the  2 pm  which  the  wave  travels  in  the  0.6  nsec 
it  takes  to  close  the  sampling  gate.  In  a system  application,  this  rolloff  can  be  eliminated  by 
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tailoring  the  response  of  the  input  transducer  and  its  matching  r stwork  to  flatten  the  overall 


response . 


D.  L.  Smythe 
B.  E.  Burke 


D.  PARALLEL-PLATE  PLASMA  ETCHING  FOR  SILICON  CCDs 

Plasma  etching  of  deposited  thin  films  of  Si^N^  and  polycrystalline  silicon  has  been  incor- 
porated into  the  fabrication  sequence  for  CCD  devices  (see  p.40  in  Ref. 2).  To  improve  process 
control,  various  strategies  have  been  investigated  to  maintain  surface  temperature  on  the  sili- 
con substrates  at  an  equilibrated  value  throughout  a particular  etching  cycle.  These  strategies 
involve  preheating  using  radiant  or  resistance  heaters,  preheating  using  an  inert-gas  plasma 
cycle,  or  constraining  the  process  to  temperatures  near  room  temperature  by  use  of  water  or 
refrigeration  cooling.  Optimum  process  control  appears  to  correspond  to  constraining  process 
temperatures  to  25°C  using  a water-cooling  coil  on  the  sample -holding  (ground)  plate. 

The  plasma-etching  system  is  similar  to  an  RF-diode  sputtering  system  (Fig.IV-8).  The 
wafers  to  be  processed  are  placed  on  top  of  a stainless-steel  ground  plate  which  comprises  one 
side  of  a capacitor  across  which  13.56  MHz  RF  power  is  applied.  The  base  plate  is  shown  in 
Fig.IV-9.  The  gas.  typically  a mixture  of  CF4  with  a low  percentage  of  02,  flows  from  an  inlet 
below  the  ground  plate  and  sweeps  radially  across  this  plate  to  a vacuum  port  in  the  center.  The 
flow  is  controlled  by  a needle  valve,  while  a mechanical  vacuum  pump  is  at  a constant  throttle. 
The  ground  plate  may  be  water-cooled  and  has  provisions  for  the  installation  of  electrical 
heaters.  Meters  for  measuring  forward  and  reflected  power  are  provided  to  assist  with  manual 


Fig.IV-8.  Parallel-plate  plasma  reactor.  At 
left  is  manifold  mounted  on  top  of  a Henry  RF 
exciter.  Vacuum  chamber  is  on  right  under  a 
tuning  network.  Below  cylindrical  vacuum 
chamber  28  cm  diam  x 14  cm  high  are  a vari- 
ety of  instruments  including  a thermocouple 
vacuum  gauge,  a capacitance  manometer,  mass 
flow  meters,  heater  variac,  gas  rotometer, 
valves,  and  regulators. 
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Fig.IV-9.  View  of  open  vacuum  chamber  with  5-cm-diam  wafer  placed 
on  23-cm-diam  base  plate.  Gas  is  exhausted  through  port  in  center  of 
stainless-steel  plate.  Gas  inlet  is  located  under  base  plate.  When 
closed,  distance  between  plates  is  8.25  cm. 


Fig.  IV-10.  Pictorial  representation  of  surface  profiler  thickness  mea- 
surements as  a function  of  position  on  a typical  plasma-etch  test  run. 
Thin  film  in  this  case  is  Si3N4,  and  negative  resist  was  utilized  earlier 
to  form  a cross -grid  etch  mask.  For  this  experiment,  ground  plate 
i was  heated  to  41  'C  using  electrical  cartridge  units  before  8-min.  etch 

cycle  commenced. 
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impedance  matching,  and  continual  adjustments  are  found  to  be  necessary  during  the  first  min- 
utes of  etching.  Typical  operating  parameters  are  a pressure  of  180  to  280  mTorr,  an  RF  power 
of  100  W forward,  a gas  flow  of  22  seem,  and  an  anode -to-cathode  spacing  of  8.25  cm.  The  sys- 
tem volume  is  8600  cm,3  and  the  ground  plate  is  23  cm  in  diameter. 

Data  were  collected  on  oxidized,  5-cm-diam  (100)  silicon  wafers  which  were  patterned  with 
negative  photoresist  on  top  of  nominal  1000-A  SijN^  or  polysilico'n  films  deposited  using  low- 
pressure  chemical-vapor-deposition  (LPCVD)  reactors.  The  photoresist  pattern  consisted  of 
lines  crossing  the  wafers  both  in  a direction  parallel  to  the  (110)  flat  and  in  a direction  perpen- 
dicular to  the  flat.  Measurements  were  performed  on  the  two  perpendicular  axes  that  cross  -the 
center  of  a wafer  using  a step  profiler.  Up  to  18  steps  were  measured  on  each  sample,  permit- 
ting determination  of  an  average  etch  rate.  Results  of  a typical  measurement  are  given  in 
Fig.IV-10.  In  addition  to  indicating  etch  rate,  these  measurements  are  useful  in  current  work 
on  determining  etch  uniformity  as  a function  of  electrical,  temperature,  gas-flow,  vacuum,  and 
geometrical  considerations. 

Plasma  etching  at  elevated  temperatures  is  attractive  because  etch  rates  increase  as  an  ex- 
ponential function  of  temperature.10  However,  this  work  indicates  that  the  higher  temperature 
process  introduces  process  control  difficulties  and  is  a factor  in  the  etch-uniformity  difficulty. 
With  no  attempt  to  control  wafer  temperature  during  etching,  a plasma-etch  run  at  100  W results 
in  a rise  in  surface  temperature  from  room  temperature  to  an  equilibrium  of  40  ° to  45  °C  in  a 
time  on  the  order  of  5 min.  Because  of  the  exponential  dependence  of  etch  rate  on  temperature,10 
there  is  no  etch-rate  control  during  this  heatup  period.  Use  of  electrical  resistance  or  radiant 
heating  on  the  base  plate  to  maintain  a constant  temperature  is  complex,  as  this  heat  input  will 
have  to  be  decreased  as  the  plasma  reaction  proceeds. 

An  alternative  is  to  try  to  preheat  the  wafers  and  ground -plate  surface  in  an  inert-gas  plasma, 
typically  nitrogen.  The  nitrogen-plasma  preheating  cycle,  however,  is  not  effective  as  a means 
of  achieving  an  equilibrated  etch  temperature.  It  is  also  more  time  consuming  than  a cooling 
process  because  the  heating  rate  is  low  and  involves  the  use  of  nitrogen  gas  which  must  be  pumped 
away  before  the  etching  gas  is  introduced,  and  thus  any  time  gained  in  etch  rate  at  higher  tem- 
peratures is  lost  by  the  length  of  the  preheating  process.  The  problem  is  graphically  illustrated 
in  Fig.IV-11,  which  indicates  the  temperature  of  the  surface  of  the  ground  plate  with  time  during 


Fig.IV-li.  Standard  process  with  N2  pre- 
heat step.  N2  plasma  is  on  for  5 min.,  dur- 
ing which  time  base -plate  surface  reaches 
51  °C.  It  takes  about  4 min.  more  to  ex- 
haust N2  from  system,  during  which  time 
base-plate  temperature  returns  to  room- 
temperature  range.  As  a result,  no  bene- 
fit is  achieved  by  N2 -plasma  preheat  step. 
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a typical  process  initiated  with  a N2 -plasma  heatup  step.  The  effect  of  the  thermal  mass  of  the 
ground  plate  is  such  that  the  net  increase  in  bulk  temperature  of  the  plate  after  5 min.  of  a \2 
plasma  is  only  4°C,  resulting  in  rapid  cooling  of  the  plate  surface  during  N2  pumpout. 

By  igniting  the  plasma  with  the  base  plate  at  22  °C  and  running  water  at  25  gph  through  the 
ground -plate  cooling  coil,  we  found  that  the  plate  surface  temperature  did  not  increase  by  more 
than  3 °C  on  a typical  run.  The  operator  needed  only  to  verify  the  starting  temperature  (20°  to 
25 °C)  and  the  final  temperature  after  the  plasma  was  turned  off.  The  water  was  left  running 
for  some  minutes  after  the  etch  run  so  that  the  temperature  would  return  to  its  original  level, 
thereby  avoiding  different  initial  temperatures  for  successive  runs.  The  average  etch  rate  of 
LPCVD  SijN^  at  22  cc/min.  and  100  W is  reduced  by  nearly  30  percent  when  the  plate  tempera- 
ture is  controlled  by  water  cooling  to  25  °C,  as  compared  with  the  etch  rate  near  the  equilibrated 
plasma-surface  temperature  of  40°  to  45 °C  which  is  achieved  after  5 min.  of  noncontrolled  op- 
eration. We  also  observed  that  the  uniformity  of  the  etching  is  better  at  lower  temperatures, 
though  temperature  is  not  the  only  contributing  factor  in  the  uniformity  problem. 

In  conclusion,  operation  of  a parallel-plate  plasma  system  with  controllable  etch  rates 
appears  to  be  optimized  by  constraining  the  temperature  of  the  base  plate  to  near  room  tempera- 
ture using  a water-cooling  coil.  Run-to-run  reproducibility  has  been  achieved  in  the  etching  of 
various  thicknesses  of  LPCVD  Si^N^  and  polysilicon,  even  when  making  the  runs  in  quick 
succession. 

D.J.  Silversmith 
R.W.  Mountain 

O.  Sotomayor-Diaz 
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V.  SURFACE-WAVE  TECHNOLOGY 


A.  PARAMETRIC  INTERACTIONS  IN  ACOUSTOELECTRIC  SURFACE 

ACOUSTIC  WAVE  (SAW)  DEVICES 

Several  types  of  devices  which  employ  nonlinear  acoustoelectric  interactions  have  been 
developed.  The  earliest  device  was  the  acoustoelectric  convolver,1  which  has  since  reached 
an  advanced  stage  of  development.  ’ More-recent  devices  — memory  correlators,  ' 3 coherent 
integrators,  and  integrating  correlators^  — have  incorporated  charge  storage  in  a convolver- 
like structure  to  effect  a variety  of  signal-processing  functions.  In  extending  this  work  on 
charge-storage  acoustoelectric  devices,  it  has  been  recognized  that  the  overall  device  perfor- 
mance results  from  a two-step  parametric  interaction.  The  first  parametric  interaction  re- 
sults in  a stored  charge  pattern,  and  the  second  interaction  utilizes  this  stored  pattern  to  pro- 
duce an  output.  The  specific  interactions  used  in  each  step  combine  to  determine  the  overall 
signal-processing  function  performed  with  the  device. 

Examination  of  the  various  two-signal  processes  for  charge  storage  and  readout  shows  that 
the  alternative  parametric  interactions  can  be  straightforwardly  interrelated.  It  has  also  been 
possible  to  extend  these  results  to  high-order  interactions  and  to  enumerate  all  possible  signal- 
processing functions  which  can  be  performed  with  this  structure.  This  report  describes  this 
functional  summary  of  parametric  interactions  in  acoustoelectric  signal  processing. 

The  structure  analyzed  here  consists  of  a piezoelectric  substrate  and  a strip  of  silicon  held 
in  close  proximity  to  the  piezoelectric  surface.  The  electric  fields  generated  by  surface  waves 
on  the  piezoelectric  interact  with  carriers  in  the  silicon.  Diodes  on  the  silicon  provide  a means 
for  integrating  and  storing  a charge  pattern  generated  by  nonlinear  interactions  in  the  silicon. 

The  back  side  of  the  silicon  is  covered  by  a uniform  conductor  or  plate.  There  exist  other 
structures  with  analogous  physical  interactions  which  can  also  be  described  within  the  frame- 
work of  this  report.  For  simplicity,  the  diode  structure  is  considered  here  as  a prototype  for 
these  other  devices. 

The  parametric  interaction  of  two  signals  which  results  in  a stored  charge  pattern  is  often 
quite  complex.  For  the  rest  of  this  report  it  will  be  useful  to  consider  a simple  model  for  the 
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charge-storage  mechanism,  which  will  allow  straightforward  examination  of  the  signal- 
processing possibilities.  In  this  model  it  is  assumed  that  the  rate  of  change  of  stored  charge 
at  any  point  on  the  silicon  is  proportional  to  the  local  product  of  the  interacting  signal  amplitudes. 

Figure  V-l(a-b)  shows  the  two  distinct  configurations  which  may  be  used  to  create  a charge 
pattern  as  the  result  of  the  application  of  two  signals;  both  have  the  form  of  a cross-correlation 
of  the  two  signals. 

The  two  possible  readout  configurations,  both  of  which  are  convolutions  with  the  stored  pat- 
tern, are  shown  in  Fig.  V-2(a-b).  Two  alternate  configurations  obtained  by  spatially  reversing 
Figs.  V-2(a)  and  (b)  cannot  be  considered  as  distinct  since  they  are  simply  convolutions  with  a 
spatially  reversed  charge  pattern,  and  this  spatial  reversal  could  be  achieved  by  interchanging 
the  original  charging  signals.  For  emphasis,  however,  the  convolution  with  the  reversed  pat- 
tern is  indicated  because  of  its  importance  as  a matched-filtering  or  correlation  process. 

Possibilities  resulting  from  the  charging  and  readout  processes  are  indicated  in  Table  V-l. 
The  two  steps  are  a correlation  of  two  signals  followed  by  a convolution  with  a third.  The 
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Fig.  V-l(a-b).  Charge  storage. 
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Fig.  V-2(a-b).  Charge  readout. 


TABLE  V-l 

SIGNAL-PROCESSING  FUNCTIONS  OBTAINABLE 

WITH  CORRELATION  OF  TWO  SIGNALS 

FOLLOWED  BY  CONVOLUTION  WITH  A THIRD  SIGNAL 

A 

B 

C 

Result 

H(t) 

Impulse 

Input 

Programmable  filter 
(memory  convolver) 

Impulse 

H(t) 

Input 

Programmable  matched  filter 
(memory  correlator) 

Pulse  train 

Impulse  train 

Impulse 

Coherent  pulse  overlay 
(coherent  integrator) 

Impulse  train 

Pulse  train 

Pulse 

Burst-waveform  processor 

Input 

Reference 

Impulse 

Correlation  receiver  with  window 
(integrating  correlator) 

Reference 

Input  + H(t) 

H(t) 

Integrating  correlator  with 
improved  readout 
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situation  of  making  the  readout  convolution  look  like  a correlation  can  be  achieved  by  several 
possible  changes  in  geometry,  but  it  is  assumed  for  this  discussion  that  it  is  achieved  by  inter- 
changing the  charging  signals. 

By  making  one  input  approximately  an  impulse,  it  is  possible  to  store  a signal  H(t)  applied 
to  the  other  input.  The  device  then  can  appear  to  be  a filter,  with  H(t)  as  its  impulse  response 
for  a subsequent  signal  C(t).  If  C(t)  also  approximates  an  impulse,  then  the  device  acts  as  an 
analog  memory  for  the  waveform  H(t).  By  interchanging  the  original  signals,  it  is  possible  to 
make  a matched  filter  for  H(t)  or  a memory  which  can  read  out  the  time  reverse  of  H(t).  In  the 
matched-filter  mode,  the  device  is  referred  to  as  a memory  correlator.^’9 

If  a pulse  train  is  one  input  and  a synchronous  impulse  train  is  the  other,  then  the  resulting 
charge  pattern  will  have  all  the  pulses  overlayed  in  the  same  location  within  the  device.  A 
device  called  a coherent  integrator6  has  been  developed  for  implementing  this  process  for  im- 
proving the  signal-to-noise  ratio  of  a repeated  signal.  An  impulse  is  used  to  read  out  the  re- 
sulting pattern.  If  the  pattern  is  not  a simple  pulse  but  has  significant  structure,  then  correla- 
tion readout  can  be  used  to  compress  the  waveform  on  readout.  The  resulting  use  of  a coherent 
integrator  is  ideally  suited  to  processing  radar  burst  waveforms. 

If  the  inputs  are  a long-duration  signal  and  a suitable  reference,  then  the  resulting  charge 
pattern  is  a finite  portion  of  the  cross-correlation  function,  and  this  can  be  read  out  to  look  for 
a correlation  spike  which  would  appear  if  the  timing  of  the  signal  and  reference  was  correct  to 
within  the  propagation  time  through  the  structure.  A device  called  an  integrating  correlator9 
has  been  developed  to  act  as  a correlation  receiver  which  searches  a large  timing  window  on 
each  integration  interval,  thus  significantly  reducing  search  time  for  acquisition  of  a spread - 
spectrum  signal. 

All  the  examples  given  so  far  have  used  an  impulse  for  at  least  one  of  the  three  applied  sig- 
nals. A more-general  capability  is  available  if  an  impulse  is  not  used.  The  last  example  of 
Table  V-l  shows  this  concept  applied  to  the  integrating  correlator.  Ralston9  has  shown  that  the 
integrating  correlator  can  be  improved  by  spatially  spreading  the  stored  correlation  function, 
by  passing  the  input  through  a dispersing  filter,  and  then  recompressing  on  readout.  But  this 
corresponds  to  replacing  the  impulse  readout  with  a cross-correlation  readout  using  the  impulse 
response  of  the  dispersing  filter. 

It  is  instructive  to  consider  the  tv  -k  diagram  representation  of  parametric  interactions. 
Proper  interpretation  of  the  tv-k  plane  as  the  domain  of  Fourier-transform  space  for  the  vari- 
ous signals  will  lead  to  correct  conclusions  concerning  these  interactions. 

Figure  V-3(a)  shows  the  tv-k  plane,  with  heavy  lines  indicating  the  allowed  dispersion  lines 
for  signals  which  occur  in  the  acoustoelectric  charge-storage  devices.  Signals  on  the  semi- 
conductor strip  must  be  spatially  uniform  and  therefore  must  fall  on  the  tv -axis.  A stored 
charge  pattern  is  at  tv  = 0 and  thus  appears  on  the  k-axis.  The  line  tv  = kv  (where  v is  the  SAW 
velocity)  represents  forward -traveling  acoustic-wave  signals,  while  tv  = -kv  represents  waves 
in  the  reverse  direction.  The  filled  and  empty  dots  represent  the  "normal"  and  "conjugate" 
components  of  a real  wave  traveling  in  the  forward  direction.  The  <v-k  diagram  can  be  used  as 
a map  of  the  two-dimensional  Fourier-transform  space  of  a general  space-time  signal.  A line 
segment  on  the  appropriate  dispersion  line  is  used  to  represent  a wideband  signal. 

These  ideas  are  now  applied  to  the  interaction  of  two  acoustic -wave  signals,  as  shown  in 
Fig.  V-3(b).  A parametric  interaction  is  represented  as  vector  addition  in  the  <v-k  plane:  the 
vector  sum  indicates  the  tv  and  k of  the  signal  generated  as  a result  of  the  interaction.  The 
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"k-shifted"  version  of  the  degenerate  case  when  the  acoustic-wave  signals  have  the  same  center 
frequency;  it  still  stores  the  correlation  of  the  two  signals  as  a charge  pattern. 

A three -signal  interaction  has  been  introduced  in  the  tv-k  discussion  to  show  compatibility 
with  the  Fourier-transform  interpretation  of  u)-k  diagrams.  The  u>-k  diagram  of  Fig.  V-3(a) 
shows  four  possible  dispersion  lines.  In  an  interaction  of  order  n,  an  output  signal  is  produced 
by  n input  signals,  and  n + 1 dispersion  lines  must  be  used.  The  number  of  distinct  interac- 
tions which  uso  only  one  dispersion  line  at  a time  is  quite  restricted.  There  are  no  fourth-  or 
higher-order  interactions  which  are  distinct  from  second-  or  third-order  possibilities  in  struc- 
tures of  this  type  for  lack  of  more  dispersion  lines. 
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Fig.  V-4(a-c).  Three-signal  interactions. 
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Figure  V-4(a-c)  summarizes  the  distinct  three-signal  interactions.  There  is  only  one 
three-signal  charge-storage  interaction,  shown  in  Fig.  V-4(a).  Since  all  three  signals  appear 
forward  in  time  inside  the  integral,  this  is  analogous  to  correlation  and  one  might  call  this 
"triple-product  correlation."  Figure  V-4(b)  shows  one  of  the  possible  readout  interactions.  It 
has  the  form  of  a convolution  of  the  two  acoustic  signals,  but  the  stored  charge  pattern  acts  as 
a stationary  weighting  function.  Others  have  referred  to  this  operation  as  "triple -product  con- 
volution." Figure  V-4(c)  shows  the  other  distinct  readout  which  has  the  form  of  a correlation 
of  the  time  signals  with  the  stored  charge,  again  as  a weighting  function.  By  analogy  with 
Fig.  V-4(b),  this  process  is  called  "triple -product  correlation." 

These  three-signal  interactions  are  interesting  for  two  reasons.  First,  if  the  third  signal 
is  a CW  pump,  then  the  interaction  can  be  a convolution  or  correlation  of  two  signals  as  in  a 
two-signal  interaction,  but  with  a shift  of  the  k of  the  stored  pattern  or  the  cv  of  the  output  sig- 
nal. These  degrees  of  freedom  might  be  of  use  in  the  design  of  practical  devices  — for  example, 
in  the  suppression  of  spurious  signals.  Second,  if  used  in  their  full  generality,  these  three- 
signal  interactions  offer  new  signal-processing  capabilities  such  as  beam  forming  with  arrays 
of  arbitrary  geometry.*** 
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Fig.  V-5.  Acoustoelectric  signal-processing  operations. 

Figure  V-5  shows  a summary  of  the  real-time  interactions,  two-signal  storage  and  readout, 
and  three-signal  storage  and  readout.  The  familiar  acoustoelectric  convolver  and  a correlation 
device1*  constitute  the  possible  real-time  interactions.  Two-signal  charge  storage  is  always  a 
correlation,  and  two-signal  readout  is  always  a convolution  process,  although  the  ability  to 
access  the  spatially  reversed  charge  pattern  can  make  this  readout  appear  as  a correlation. 
Three-signal  charge  storage  has  the  form  of  a triple -product  correlation  or  a k-shifted  correla- 
tion of  two  signals.  Three-signal  readout  is  either  triple -product  convolution  or  triple -product 
correlation,  but  can  also  be  used  to  perform  u -shifted  convolution  or  correlation  of  two  signals. 

J.  H.  Cafarella 
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B.  RESONANT  INTERACTIONS  OF  SURFACE  AND  BULK  ACOUSTIC 

WAVES  IN  GRATINGS 

In  transmitting  a surface  wave  through  normal-incidence  gratings  with  a period  approxi- 
mately equal  to  the  acoustic  wavelength,  one  observes  a series  of  sharp  stop  bands,  as  shown 

o o 

on  Fig.  V-6  for  gratings  of  two  different  depths,  1920  and  3640  A.  For  a 4800-A-deep  grating, 
the  stop  bands  are  deeper  and  sharper.  Expanded  plots  of  the  transmission  at  one  such  stop 
band  are  shown  in  Fig.  V-7. 

The  stop  band  occurs  when  the  surface  wave  couples  to  a plate  mode  of  the  substrate.  The 
modes  to  which  the  surface  wave  couples  consist  essentially  of  bulk  waves  which  are  multiply 
reflected  at  near-vertical  incidence  between  the  top  and  bottom  surfaces.  The  condition  for 
coupling  can  be  expressed  as 

Kp  = Kg-Ks  (V-l 

where  K is  the  projection  along  the  surface  of  the  wave  vector  of  the  plate  mode,  K = 2ir/d 
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Fig.  V-6.  Transmission  through  a 200-groove  Fig.  V-7.  Same  as  Fig.  V-6  except  expanded 

grating  at  normal  incidence.  Center  frequency  . f 4800-A-deen  ffratimr  Center  fre- 

is  170  MHz.  Period  of  grooves  is  20.32  am.  ana  tor  a 4BUU  A oeep  grating.  center  Ire 

6 ^ quency  is  194.2  MHz. 


where  d is  the  grating  period,  and  K is  the  magnitude  of  the  wave  vector  of  the  surface  wave 

s 

in  the  grating.  This  is  simila.'  to  the  problem  of  coupling  of  interdigital  transducers  to  plate 
modest 

A simplified  model  of  the  coupling  can  be  derived  by  first  examining  the  angle  of  propaga- 
tion of  a scattered  bulk  wave,  and  then  specifying  the  conditions  for  efficient  coupling  into  plate 
modes  characterized  by  a bulk  wave  propagating  at  this  angle.  The  exact  solution  for  a plate 
mode  involves  both  shear  and  longitudinal  bulk  waves.  In  order  to  derive  a simplified  model  of 
the  resonance  phenomenon,  the  approximation  is  made  that  only  a single  bulk  wave  is  scattered 
and  that  there  is  no  longitudinal-to-shear-wave  coupling  upon  reflection.  With  these  approxi- 
mations Eq.  (V-i)  can  be  rewritten  as 


|Kb|  cose  = Kg  - Ks 


(V-2) 


where  KB  is  the  wave  vector  of  the  only  bulk  wave  in  the  model,  and  e is  the  angle  relative  to 

the  surface  at  which  the  bulk  wave  is  scattered  and  reflected  (see  Fig.  V-8).  When  K = K or, 

g s 

equivalently,  when  the  wavelength  of  the  surface  wave  is  approximately  equal  to  the  grating 
period,  scattering  occurs  at  8 = 90°. 

Equation  (V-2)  may  be  rewritten  as 


cos  8 = 


( V-3) 


where  v_  and  v are  the  velocities  of  the  bulk  and  surface  waves,  and  f = v /d.  Note  that  f 

a o S O 

is  twice  the  frequency  of  the  stop  band  in  which  coupling  of  a surface  wave  into  a reflected  sur- 
face wave  occurs.  For  f > f the  bulk  waves  are  scattered  in  the  reverse  direction  (8  < 90°), 
while  for  f < f the  scattering  is  in  the  forward  direction. 


TT-t-istrti 
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Fig.  V-8.  Model  of  sharp  transmission 
stop  bands,  r ^ and  t3  are  bulk -wave 
transit  times,  and  Tj  is  surface-wave 
transit  time  from  conversion  points, 
n is  an  integer. 


A stop  band  is  produced  when  the  bulk  wave  bounces  off  the  bottom,  travels  back  up  to  the 
grating,  and  reconverts  to  a surface  wave  which  is  out  of  phase  with  the  original  surface  wave 
traveling  in  the  grating.  This  is  illustrated  in  Fig.  V-8.  When  many  reflections  occur,  a sharp 
stop  band  results.  Therefore,  these  sharp  stop  bands  occur  only  when  8 ~ 90°,  i.e.,  when 


Several  experimental  observations  indicate  the  existence  and  nature  of  the  surface-wave- 
to-plate-mode  coupling  in  normal-incidence  gratings.  When  a surface  wave  is  launched  on  the 
top  (grating)  surface  of  a substrate,  the  spatial  period  of  a plate  mode  may  be  measured  by 
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surface -probing  the  field  pattern  on  the  bottom  surface.  This  was  done  for  a series  of  frequen- 
cies corresponding  to  the  centers  of  several  stop  bands.  In  each  case,  the  measured  periodicity 
was  accurately  predicted  by  Eq.  (V-l). 

In  a second  series  of  experiments,  a CW  tone  was  applied  to  the  input  transducer  and  the 
result  of  multiple  bounces  of  the  scattered  waves  was  observed  vs  time  in  the  output  transducer, 
as  shown  in  Fig.  V-9.  The  top  left  photo  shows  the  received  RF  signal.  The  next  three  photos 
show  the  phase-detected  response  for  frequencies  at  the  middle  of  a stop  band,  between  two  stop 
bands,  and  at  the  middle  of  the  next  stop  band,  respectively.  Note  that  the  multiple  echoes  are 
out  of  phase  with  the  direct  transmission  at  the  centers  of  the  stop  bands,  whereas  the  echoes 
alternate  in  phase  at  a frequency  between  stop  bands. 

As  a third  check  on  the  existence  of  plate-mode  coupling,  the  back  of  a substrate  was 
scored  with  a diamond  saw;  the  multiple  stop  bands  disappeared,  as  anticipated. 


Fig.  V-9.  Transmission  of  a 200-nsec  tone  burst  through  grating.  Upper  left  is  direct  RF 
output.  At  top  right,  it  is  shown  mixed  with  CW  to  yield  phase  information  at  a 185.99- MHz 
stop-band  frequency.  Left  bottom  trace  is  between  two  stop-band  frequencies  at  188.08  MHz 
showing  that  alternate  bounces  of  bulk  wave  yield  out-of-phase  signals.  Bottom  right  trace 
is  at  next  stop-band  frequency  (190.02  MHz). 
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Using  the  condition  for  resonance  shown  in  Fig.  V-8,  one  can  make  an  approximate  calcu- 
lation of  the  stop-band  frequencies.  The  timing  condition  of  Fig.  V-8  combined  with  Eq.  (V- 3) 
yields  the  relation 


(n  + l>  7 

2bVvB 


(V-4) 


where  b is  the  thickness  of  the  substrate.  The  quantity  ZbfQ  /v^  is  the  number  of  bulk-wave 
periods  in  the  round  trip  from  the  top  to  the  bottom  of  the  crystal  at  normal  incidence  at  fre- 
quency fQ.  Equation  (V-4)  may  be  solved  for  the  stop-band  frequencies  f corresponding  to  suc- 
cessive values  of  n.  The  parameters  inserted  into  Eq.  (V-4)  to  calculate  the  values  of  f for 
the  experiments  on  Y-Z  LiNbO-j  were  vg  = 3485  m/sec,  b = 1.27  mm,  and  fQ  = 171.51  MHz.  In 
addition,  the  dominant  bulk  wave  was  assumed  to  be  a shear  wave  propagating  at  a constant 
velocity  of  4444  m/sec  (Ref.  13). 


5 180 ! 


llt-8-15?IOl 


/ 


A 

/• 


/ 

/• 

/• 


4 

/ 

i 

o • 


• / 

</ 


• MEASURED 
-o-  CALCULATED 


10 

NUMBER 


Fig.  V-10.  Calculated  and  measured  stop-band  frequencies  plotted 
in  order  vs  "number."  Measured  stop  bands  are  those  obtained  from 
a display  such  as  shown  in  Fig.  V-6.  Stop  bands  plotted  are  those 
which  are  most  prominent  and  appear  to  belong  to  a series.  When 
a stop  band  appears  split,  several  frequencies  are  plotted  at  one 
number. 


The  theoretical  stop-band  frequencies  and  the  frequencies  of  the  most -prominent  experi- 
mental stop  bands  are  plotted  in  Fig.  V-10.  The  correspondence  is  seen  to  be  good  for  f > fQ. 
(The  curves  can  be  displaced  laterally  since  the  "number"  assignment  is  arbitrary.)  The  dis- 
agreement for  f < fQ  may  be  due  to  the  approximations  used.  The  deviation  of  the  simple  model 
is  in  the  right  direction,  since  the  inclusion  of  longitudinal  waves  in  the  exact  plate-mode  calcu- 
lation would  yield  higher  coupled-wave  velocities  and  a lower  slope  for  the  theoretical  curve. 
However,  the  approximate  quantitative  agreement  indicates  that  the  general  interpretation  of 
the  phenomenon  is  correct. 


J.  Melngailis 
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